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Introduction 


As  an  alternative  to  electric  power  generation  by  non-renewable  energy 
sources  such  as  gas,  oil  and  coal,  geothermal  brines  can  be  used  as  heat 
sources  in  power  cycles.^  In  one  design,  the  dry  steam  flashed  from  the 
brine  is  expanded  in  a turbine.  Since  geothermal  temperatures  are  much 
lower  than  those  used  in  conventional  steam  power  generation,  the  use  of  a 
working  fluid  with  a critical  temperature  much  lower  than  that  of  steam  has 
been  advocated.  Pure  isobutane  was  the  first  choice  of  working  fluid.  By 
using  a mixture  of  isobutane  and  isopentane,  however,  two  advantages  are 
obtained:  the  location  of  the  critical  point  can  be  optimized  by  modifying 

the  composition;  and  the  performance  of  the  heat  exchanger  is  improved 
because  the  mixture,  contrary  to  the  pure  fluid,  undergoes  the  isobaric 
phase  transition  from  liquid  density  to  subcritical  density  over  a range 
of  temperatures,  rather  than  abruptly.  For  these  reasons,  the  demonstration 
project  which  is  being  constructed  under  sponsorship  of  the  Department  of 
Energy  will  use  a 90  mol  % isobutane  - 10  mol  % isopentane  mixture  as 
a working  fluid. 

In  the  design  of  the  power  cycle,  it  is  necessary  to  have  information 
on  the  thermodynamic  and  transport  properties  of  the  working  fluid.  For 
the  past  six  years,  groups  at  the  National  Bureau  of  Standards  have  worked 
at  obtaining  the  data  base,  validating  existing  data  and  models  and  developing 
correlative  and  predictive  methods  for  establishing  these  property  values. 

This  report  contains  the  data  obtained  under  this  contract,  their  evaluation, 
and  the  predictive  methods,  previously  existing  ones  as  well  as  those  developed 
for  this  project.  The  report  consists  of  three  parts,  two  referring  to  the 

thermodynamic,  the  third  to  the  transport  properties.  The  approaches  used  in 

i 

the  two  cases  are  somewhat  different.  In  the  case  of  the  thermodynamic  properties, 
no  predictive  methods  were  available  even  for  pure  isobutane,  let  alone  for 
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the  mixtures.  Thus,  in  the  early  stages  of  the  work,  the  vapor  pressure 

and  several  PV  isotherms  of  supercritical  and  liquid  isobutane  were 
2 3 4 

measured.  ’ ’ This  enabled  NBS  scientists  to  select  the  reliable  data 

from  the  body  of  data  available  in  the  literature, ^ On  the  basis  of 

these  validated  data,  a thermodynamic  surface  was  constructed.  This 

surface  was  not  accurate  in  the  critical  region.  On  the  basis  of  the 

Wilson-Kadanoff  renormalization  group  theory,  a separate  nonclassical 

(scaled)  surface  was  constructed  for  the  critical  region.^ 

For  the  isobutane-isopentane  mixtures,  the  first  step  was  to  obtain 

(P,x,T)  data  for  the  critical  line.  In  the  process,  the  critical  point 

of  pure  isopentane,  and  vapor  pressure  data  between  370  K and  the  critical 

point  (408  K)  were  measured.  In  addition,  the  pressures,  densities  and 

compositions  of  coexisting  phases  of  mixtures  of  approximately  50  mol  % 

isobutane,  50  mol  % isopentane  were  measured  in  the  range  of  290-325  K. 

Within  the  framework  of  another  contract,  a global  thermodynamic  surface 

was  developed  on  the  basis  of  the  principle  of  generalized  corresponding 

states.  In  part  II  of  this  report,  the  data  for  pure  isopentane  and  for 

the  critical  line  are  used  to  develop  a nonclassical  (scaled)  thermodynamic 

surface  for  the  mixtures.  This  surface  was  intended  to  be  valid  within 

a range  of  approximately  3 K below  to  30  K above  the  critical  line.  We 

will  compare  the  scaled  surface  for  the  mixture  with  the  global  surface 

8 9 

reported  elsewhere.  ’ Unfortunately,  the  agreement  is  not  good.  Our 
conclusion  is  that  the  scaled  surface  needs  further  basic  work.  We 
have,  however,  developed  an  accurate  scaled  model  for  representing 
coexisting  phases  of  isobutane-isopentane  mixtures. 

The  work  on  the  transport  properties  has  been  developed  from  a 
different  perspective.  Predictive  techniques  for  the  transport  properties 
of  hydrocarbon  fluids  were  available  at  the  onset  of  this  project  in  the 
form  of  the  program  TRAPP  developed  by  NBS  scientists'^’ on  the  basis  of 
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generalized  corresponding  states.  The  experimental  program  on  transport 

properties  was  therefore  directed  at  validation  of  the  TRAPP  program  for 

pure  fluids  and  mixtures,  over  as  large  a range  of  reduced  temperatures  and 

pressures  as  experimentally  feasible.  Since  the  available  viscometer  has  an 

12-17 

upper  temperature  limit  of  320  K,  initially  cryogenic  fluids  and  some 

18“20 

of  their  mixtures  such  as  nitrogen,  methane  and  ethane  were  studied. 

This  work  led  to  the  conclusion  that  the  program  TRAPP  predicts  viscosities 

correctly  on  the  + 5%  level  for  densities  smaller  than  2p^. 

In  the  second  half  of  the  work,  however,  experimental  studies  were 

21 

done  of  the  two  isomers  normal  butane  and  isobutane.  It  was  then  found  that 
there  is  a rather  substantial  isomeric  effect  in  the  viscosities  of  the  liquid 
that  is  not  well  predicted  by  the  program  TRAPP.  An  alternative  formulation 
of  the  viscosity  of  liquid  isobutane  is  proposed.  This  study  of  the 
transport  properties  of  hydrocarbons  has  enabled  us  to  identify  the 
ranges  of  experimental  parameters  in  which  the  results  of  the  program 
TRAPP  can  be  used  with  confidence. 

In  the  first  two  parts  of  this  report,  we  will  present  the  experimental 
thermodynamic  data  obtained  under  this  contract,  and  give  the  results  of  the 
correlations  developed  with  estimates  of  their  reliability.  In  the  third 
part,  we  will  do  similarly  for  the  transport  properties,  and  compare  with 
the  prediction  of  the  program  TRAPP.  In  Appendix  A,  we  give  a complete 
listing  of  the  publications  that  have  resulted  from  this  project. 
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Part  I 


Thermodynamic  Properties  of  Isobutane  and  Isobutane-Isopentane  Mixtures' 

1.  Experimental  data  obtained  under  this  contract 

1.1  Summary  of  Burnett  PVT  data  for  isobutane 

In  a high-quality  Burnett  PVT  apparatus,  expansions  were  performed  along 

two  isotherms  in  the  supercritical  fluid  and  along  two  isotherms  in  the 

subcritical  vapor  for  a 99.98%  pure  sample  of  isobutane.  In  addition,  the 

vapor  pressure  of  isobutane  was  measured  from  298.15  to  398.15  K.  All  results 

2 3 

have  been  published,  with  an  assessment  of  the  reliability  of  the  data.  ’ The 
supercritical  PVT  data  are  listed  in  Table  II,  the  virial  coefficients  at 
the  four  temperatures  in  Table  12,  and  the  vapor  pressure  in  Table  13.  The 
authors  estimate  the  inaccuracy  of  pressures  above  0.5  MPa  as  0.01%.  The 
inaccuracy  of  densities  and  compressibility  factors  is  at  most  an  order  of 
magnitude  larger. 

1.2  Critical-point  and  critical-line  data  on  isobutane,  isopentane  and 
their  mixtures 

1.2.1  Visual  Cell 

For  observation  of  the  meniscus  disappearance  in  samples  of  isobutane, 

isopentane  and  their  mixtures,  a simple  optical  cell  was  constructed  (Fig.  II). 

It  consists  of  a stainless-steel-316  ring,  of  inner  diameter  2.54  cm,  thickness 

1.27  cm,  sandwiched  between  two  circular  sapphire  windows  (of  diameter  5 cm, 

thickness  3 mm)  by  means  of  an  aluminum  clamp.  Tin  foil  of  0.25  mm  thickness 

3 

serves  as  a gasket.  This  way,  a disk-shaped  cell  of  nominally  6.5  cm  volume 

is  formed.  The  dimensions  of  the  ring  were  calibrated  at  20  °C  by  the  NBS 

Metrology  Division.  From  the  calibration,  the  volume  of  the  cell  at  20  °C 

3 

was  calculated  to  be  6.4484  cm  , with  an  uncertainty  of  0.1%;  the  thickness 


of  the  tin  foil  gasket  has  been  taken  into  account.  The  cell  contains  a 

stirrer,  a soft-magnetic-steel  rod  which  is  moved  around  by  means  of  a 

3 

hand-operated  external  magnet.  The  volume  of  the  stirrer  is  126  mm  . One 

end  of  a capillary  of  1.6  mm  outer  diameter,  0.45  mm  inner  diameter  and 

8 cm  length  is  welded  into  the  ring.  The  other  end  is  sealed  off  by  means  of 

3 

a miniature  pressure  valve.  The  combined  volume  in  line  and  valve  is  22  mm  , 

as  calculated  from  the  dimensions.  The  coefficients  of  thermal  expansion  of 

the  steel  ring,  which  determines  the  cell  volume  at  temperatures  other  than 

— fi 

20  °C,  is  assumed  to  be  16.5  x 10  per  K in  the  range  of  20-200  °C.  This 
conclusion  was  reached  on  the  basis  of  thermal  expansion  data  for  stainless 
steel  304  in  Ref.  22,  and  the  observation,  in  Ref.  23,  that  the  stainless 
steels  304  and  316  have  identical  expansion  coefficients  below  300  K. 

The  assembled  cell  and  valve  have  a mass  of  just  over  400  g.  The 
sample  mass,  typically  of  the  order  of  1.5  g,  is  determined  by  weighings 
of  the  empty  and  the  filled  cell  on  a 1-kg  balance  with  resolution  of  0.5  mg. 

The  cell  is  connected  to  a manifold  that  contains  a miniature  Sensotech* 
strain-gage  pressure  transducer  and  another  valve  (Fig.  I 1).  Sample  gas 
vented  from  the  cell  is  captured  above  water  in  an  inverted  graduated 
cylinder.  The  pressure  transducer  was  claimed  by  the  manufacturer  to  be 
precise  to  0.1%  in  a temperature  range  up  to  200  °C.  The  arrangement  of 
Fig.  1 permits  in-situ  calibration  of  the  pressure  transducer  with  respect 
to  a dead-weight  gage.  We  found  the  transducer  performance  satisfactory  at 
room  temperature;  pressure  measurements  at  temperatures  above  130  °C,  however, 
were  hampered  by  strong  drifts  of  the  transducer,  often  as  much  as  10  kPa 
(0.1  bar)  per  day;  these  measurements  are  therefore  not  as  reliable  as  we  had 
hoped  them  to  be  on  the  basis  of  the  manufacturer’s  claim.  In  the  measurements 
on  pure  isopentane,  however,  we  mounted  the  pressure  transducer  outside  the 
bath,  so  that  these  pressures  have  better  reliability  than  those  of  the  mixtures. 
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The  measurement  procedure  is  to  weigh  in  the  required  amount  of 
isopentane,  and  then  put  in  the  matching  amount  of  isobutane  by  means  of  a 
volume  pump,  so  as  to  obtain  a mixture  of  the  desired  composition  at  an 
overall  fill  density  5-10%  above  the  estimated  critical  value.  The  amounts 
of  isobutane  and  isopentane  are  obtained  by  weighing.  The  cell  is  then 
hooked  up  to  the  manifold  and  placed  inside  a commercial  "visibility  bath" 
filled  with  a low-vapor-pressure  silicone  oil;  the  temperature  of  the  bath 
is  controlled  to  a few  hundredths  K.  The  mixture  is  heated  until  the  meniscus 
disappears  from  the  cell  at  the  top.  The  cell  is  then  heated  several  degrees 
while  the  line  and  pressure  transducer  sections  are  evacuated.  The  fluid  is 
expanded  into  the  pressure  transducer  section.  Care  is  taken  that  the 
mixture  does  not  phase-separate  in  the  cell.  A full  expansion  corresponds 
with  about  5%  decrease  in  density.  It  is  not  necessary  to  do  a complete 
expansion,  although  this  was  done  most  of  the  time.  With  the  cell  valve  open, 
the  fluid  is  now  cooled  in  steps,  and  the  pressure  and  temperature  measured, 
until  a dark  brown  color  in  the  cell  betrays  the  closeness  of  a phase 
separation.  The  cell  valve  is  then  closed  and  the  temperature  lowered  until 
the  meniscus  appears.  The  meniscus  level,  at  first  appearance,  and  after 
vigorous  stirring,  is  read  with  a cathetometer . It  will  generally  be  in 
the  upper  part  of  the  cell.  The  cell  is  then  heated  until  well  into  the 
one-phase  region,  while  the  fluid  in  the  transducer  is  expanded  into  the 

3 

inverted  graduated  cylinder  (approximately  20  cm  at  room  temperature) . 

After  the  vent  valve  is  closed,  the  sequence  of  pressure  measurements  and 
meniscus  level  determination  is  repeated.  The  critical  point  temperature  and 
density/  are  obtained  by  interpolating  in  the  meniscus  level  measurements  for 
meniscus  disappearance  at  mid-level  in  the  cell  (after  allowances  for  the 
capillary,  valve  and  stirrer  volumes  are  made).  The  critical  pressure  is 
obtained  by  extrapolating  the  isochoric  pressures  to  the  transition  points, 
and  interpolating  between  isochores.  In  several  instances  we  stopped  the 
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sequence  when  the  meniscus  disappeared  at  mid-level,  and  immediately,  in  situ, 
calibrated  the  transducer  with  respect  to  a calibrated  dead-weight  gage. 

This  alleviated  somewhat  the  problems  that  were  caused  by  the  strong  drift 
of  the  pressure  transducer.  The  barometric  pressure  is  read  on  an  aneroid 
gage  with  a resolution  of  25  ym  of  mercury. 

After  all  measurements  are  completed,  the  cell  is  removed,  cleaned  and 
weighed,  after  which  the  remaining  sample  is  blown  off  and  the  empty  cell 
weighed  again.  The  amounts  released  from  the  cell  during  the  measurements 
are  calculated  from  the  gasometry  data,  corrected  for  gas  nonideality.  Proof 
of  a successful  series  of  measurements  is  that  the  gasometry  data  add  up  to 
the  weight  loss  observed.  We  have  rejected  those  runs  in  which  the  mismatch 
between  the  gasometry  and  weight  data  exceeded  15  mg.  (For  those  runs  we 
retained,  the  agreement  was  usually  an  order  better  than  that.)  Therefore,  no 
uncertainty  larger  than  1%  in  the  reported  density  results. 

The  data  obtained  with  this  apparatus  are  summarized  in  Table  I 4 (isochoric 
PpT  data) , Table  I 5 (critical-line  data) , Table  I 6 (vapor  pressure  of  pure 
isopentane)  and  Table  I 7 (dew-bubble  data).  We  estimate  the  reliability  of 
the  data  to  be  on  the  level  of  20  mK  in  temperature,  4 kPa  in  pressure,  1%  in 
density  and  0.001  in  composition. 

The  isopentane  used  in  this  work  was  Phillips*  research  grade  of  99.99+  mol 
percent  claimed  purity,  the  isobutane  was  Phillips*  research  grade,  claimed  to 
be  at  least  99.9%  pure.  The  isobutane  liquid  was  analyzed  by  chromatography 
after  the  measurements  of  the  isobutane  vapor  pressure  and  found  to  contain 
30  ppm  of  nitrogen  and  300  ppm  of  propane.  The  n-butane  content  could  not 
be  determined. 

The  isopentane  was  kept  at  atmospheric  pressure  and  must  therefore 
have  contained  some  dissolved  air.  When  introducing  it  as  the  first  component 
into  the  visual  cell,  we  always  overfilled  and  blew  off  some  vapor  by  mild 
heating  of  the  cell,  thus  eliminating  most  of  the  volatiles. 
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1.3  VLE  data 


A schematic  of  the  VLE  cell  and  manifold  is  given  in  Fig.  I 2.  The 

pressure  cell  is  a sapphire  tube;  stainless-steel  plugs  inserted  in  the  tube 

24 

at  both  ends  seal  it  by  means  of  0-rings.  Samples  are  introduced  from  a 

thermal  compressor  which  is  removable  and  is  filled  from  a gas  buret;  the 

amounts  of  the  components  of  a mixture  are  determined  by  weighing.  The 

25 

compositions  thus  obtained  are  reliable  to  + 0.0002  in  mole  fraction. 

Mercury  fills  a calibrated  injector  pump,  the  capillaries  and  part  of  the 

sample  cell  (Fig.  12).  It  is  used  to  displace  the  sample  from  the  thermal 

compressor  into  the  cell,  to  serve  as  a pressure-transmitting  medium  and  to 

fill  all  noxious  volumes,  so  that  the  mixture  sample  is  strictly  confined  to 

the  sapphire  tube.  The  bore  of  the  sapphire  tube  was  calibrated  over  its 

length  by  means  of  triple-distilled  mercury;  that  of  the  injector  pump  was 

similarly  calibrated  by  means  of  decane,  the  density  of  which  was  determined 

3 

by  pycnometry.  The  pump  is  calibrated  to  + 2 mm  . The  pressure  is  measured 
by  means  of  a Validyne*  pressure  transducer,  kept  at  ambient  temperature  and 
calibrated  with  respect  to  a dead-weight  tester.  The  transducer  is  filled 
with  mercury  in  order  to  eliminate  noxious  volumes.  Pressures  are  measured 
to  0.1%  and  head  corrections  are  made.  The  barometric  pressure  is  read  on 
an  aneroid  gage  with  a resolution  of  25  pm  of  mercury.  The  system  is 
immersed  in  a stirred  water  bath,  controlled  to  better  than  1 mK.  The 
temperature  is  read  on  a quartz  thermometer  calibrated  with  respect  to 
a standard  platinum  resistance  thermometer.  The  temperature  resolution  is 
0.3  mK,  the  accuracy  1 mK. 

VLE  measurements  are  obtained  by  measuring  the  meniscus  levels  of 
mercury  and  the  vapor-liquid  interface  of  a predetermined  amount  of  mixture 
prepared  gravimetrically , at  measured  values  of  pressure  and  temperature. 
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At  fixed  P and  T,  the  compositions  and  densities  of  coexisting  phases  are 

nonvariants.  By  introducing  at  least  two  samples  of  different  known 

compositions,  the  densities  of  the  coexisting  phases  can  be  calculated  from 

26 

the  observed  phase  volumes  and  the  known  total  mass.  In  this  way,  the  liquid 
densities  listed  in  Table  I 8 were  obtained.  In  this  temperature  range,  the 
vapor  densities  are  too  low  to  be  obtained  with  any  accuracy  by  this  method. 

The  VLB  data  have  been  proven  highly  useful  in  the  construction  of  a 
thermodynamic  surface,  which  is  under  way  under  a different  contract. 

2.  Thermodynamic  surfaces  for  Isobutane  and  for  isobutane- isopentane  mixtures 

2.1  Global  thermodynamic  surface  for  isobutane 

On  the  basis  of  validated  PVT  data  for  isobutane,  a thermodynamic  surface 

was  constructed  that  spans  the  range  of  250-600  K at  pressures  from  0.1  to  40  MPa. 

The  surface,  and  a comparison  with  data,  are  described  in  detail  in  two  publications 

2 3 

by  Waxman  and  Gallagher.  * At  given  pressure  and  temperature,  the  surface  will 

give  densities  that  are,  in  general,  reliable  to  approximately  0.1%.  The 

reliability  deteriorates  near  the  critical  point.  The  surface  was  not 

constrained  to  the  physical  critical  point  because  such  constraint  usually 

leads  to  a deterioration  in  other  areas,  such  as  in  the  supercritical  region. 

In  a small  region  around  the  critical  point,  defined  by 

0.985  < T /T  < 1.015 
c 

and  0.7  < p^/p  < 1.3 

the  surface  should  not  be  used.  In  this  report,  we  will  give  the  equations 
defining  the  surface  in  a dimensionless  form  (Refs.  2 and  3 are  dimensional). 

The  Helmholtz  function  in  this  form  can  be  expressed  as 

A - Aq(T)  + A^(V,T)  + A^CV,!)  + A^CV,^),  1(1) 

where  A^  is  the  ideal  gas  contribution,  written  in  the  form: 
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1(2) 


A 

o 


^^00  "'''^01^^ 


_ 8 
£nT  + 2 

i=2 


A 


+ AQgT^nCe^o'^^-l)  ; 


A,  V and  T are  the  dimensionless  Helmholtz  free  energy,  volume  and  temperature, 

to  be  defined  shortly.  A^^  - A^„  and  x are  dimensionless  constants  listed 

00  09  o 

in  Table  19.  A^  is  the  "hard  sphere"  contribution,  of  the  form 


Jin 


+ 


4y 


V(l-y)  2(l-y) 


1(3) 


y is  given  by 

y = b p/4  = b/4V,  1(4) 

where  b is  an  equivalent  hard-sphere  volume,  itself  a function  of  temperature. 
The  representation  of  y as  a function  of  temperature  is 


V 


Yo  y- 


yo  + y^  £nT  + + 


1(5) 


The  dimensionless  constants  y^  - y^  are  listed  in  Table  I 9, 
A_  is  the  "second  virial"  contribution,  of  the  form 


A = ^ 
^2  - 
V 


^20 


21 


- -3 

T T 


A A 

23  , 24 

T^  plO 


1(6) 


the  dimensionless  constants  A^q  - A^^  are  given  in  Table  I 9.  Finally,  the  free 
energy  contains  a "residual"  term  A^.  It  is  obtained  by  fitting  a function  of 
the  form  ^ ' 


2 2 B. .z 

• • 

1 3 


i+1 


j 


with  z 


1 

1-e  o 


1(7) 


to  the  differences  between  experimental  data  and  calculated  values  based  on  the 

sum  of  A , A, , and  A„ . The  dimensionless  constants  B. . and  z are  given  in 
o’  1’  2 ij  o ^ 

Table  I 9. 
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The  dimensionless  quantities  are  defined  as  follows: 

T = T/T*,  1(8) 

V - V/V*, 

P - P/P*, 

A ==  A/A**, 

^ = S/S**, 

where  T*,  V*,  P*  are  dimensioned  reference  values  listed  in  Table  I 9,  and 

A**,  S**  dimensioned  reduction  factors  for  free  energy  and  entropy,  resp., 
derived  from  V*,  T*  and  P*.  Also  required  is  the  reduced  gas  constant 

R**  = R/S**.  1(9) 

The  other  thermodynamic  functions  can  be  obtained  from  combinations  of  the 
derivatives  of  A with  respect  to  V and  T in  the  usual  way,  for  instance  the 
pressure 

P = - > 1(10) 

8V  T 

the  entropy 

S = “ (~)  „ 1(11) 

8T  V 

and  the  Gibbs  function 

G = A + PV  . 1(12) 


2.2  Global  thermodynamic  surface  for  isobutane-isopentane  mixtures 

A global  thermodynamic  surface  for  isobutane-isopentane  mixtures, 

based  on  the  global  surface  for  isobutane,  described  in  Sec.  2.1,  as  a 

reference,  and  on  a generalized  principle  of  corresponding  states,  is 

9 

being  developed  under  a different  contract.  The  principal  equations  and 
preliminary  results  are  summarized  here. 
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I 

1 


h°  = 

5 4 


1(13) 


The  first  step  in  the  process  is  to  construct  a thermodynamic  surface 

for  isopentane  based  on  that  of  isobutane  as  a reference.  We  define  the  ratios  jj 

f°  = "x  ^ 

-L  5 ' 4 ’ 


e°  = P ^/P  ^ 

where  the  subscript  5 refers  to  isopentane,  the  subscript  4 to  isobutane,  and 

the  superscript  c indicates  a critical-point  value,  while  the  symbol  V indicates 

a molar  volume.  The  properties  of  isopentane  are  made  dimensionless  by  means 

of  the  same  reference  constants  as  those  of  isobutane  (Table  19).  By  our 

" c — c — • c 

choice  of  reference  constants,  T^  , P^  and  happen  to  be  unity. 

The  principle  of  corresponding  states  postulates  for  the  compressibility 
factor  Z = PV/RT: 

(V3,  T^)  = Z^  (V^  = T^  = 


cf 


and  for  the  configurational  Helmholtz  free  energy  A 


(V^,  T^)  = f°  A^^^  (V^  = V^/h°,  T^  = T^/f°)  - R**T^  Inh' 


1(14) 


for  all  values  of  V,  T.  It  implies  that 

f°  = h°  g‘ 


1(15) 


Since  no  two  substances  obey  the  principle  of  corresponding  states 

27  28  29 

accurately,  Rowlinson  et  al.  and  Leland  et  al  ’ proposed  to  generalize 

the  principle  by  the  use  of  two  functions,  9(V,T)  and  cj)(V,T)  that  ensure  an 

cf 

exact  mapping  of  Z and  A from  one  substance  to  the  other. 
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We  define  these  functions  by  means  of 

f = e(?^,f^) 

and 

g = . 1(16) 

in  such  a way  that 

ZjCVj,  Tj)  = Z^(V^  = Vj/h,  = Tj/f)  , 

(Vj,  Ij)  = f = Vj/h,  = Ij/f)  -R**Tj  Inh  , 1(17) 

with  h given  by 

h = f/g  = (Z^^/Z^^)  4>(V^J^)  . 1(18) 

We  have  used  the  experimental  data  available  for  the  vapor  pressure  of 

isopentane  (Table  I 6)  and  for  coexisting  vapor  and  liquid  densities 

30 

determined  by  Young  at  the  turn  of  the  century  to  calculate  point  values 
of  and  (p  (V^,T^) . We  have  represented  9 and  as  simple  functions 

of  density  and  temperature: 

e = a^  + ^^(p^  - 1)  + ^^(T^  - 1) 

and 

0 = b^  + - 1)  + b3(T^  - 1)  . 1(19) 

The  values  of  the  constants  a^  - b^  are  given  in  Table  I 10.  The  constants 
a^,  b^  have  been  set  equal  to  unity  because  their  fitted  values  did  not 
depart  significantly  from  unity. 

The  thermodynamic  surface  of  isopentane  is  completely  defined  by: 

the  reference  Helmholtz  function  for  isobutane  (Table  I 9) ; the  ideal-gas 

Helmholtz  function  for  isopentane,  obtained  by  interpolation  in  the  table 
31 

of  Scott  ; the  shape  factors  9 and  <})  characterized  by  the  constants  listed 
in  Table  I 10;  and  the  critical  constants  of  isobutane  and  isopentane  listed 
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in  Table  I 11.  Since  the  critical  temperature  of  the  global  isobutane  surface 


is  about  2 K above  the  physical  critical  point,  a similar  offset  will  occur  in 
the  isopentane  surface. 

The  next  step  is  the  development  of  a thermodynamic  surface  for  the 
mixtures  of  isobutane  and  isopentane.  This  is  achieved  by  again  invoking  the 
principle  of  corresponding  states.  For  a mixture  of  composition  x,  molar 
volume  V and  temperature  T we  obtain  the  molar  configurational  Helmholtz 
free  energy  from  that  of  isobutane  by  the  relations 


and 


I 

X 


1(20) 


Here  h , f are  defined  by 

X X ■' 


’'x  ' [1  + x{e(V^,  T^)  - 1}] 


and 


- 1}]. 


1(21) 


Left  to  be  defined  are  the  pseudocritical  constants  V_ 


c 


T 


c 


in  1(21).  We 


X 


X 


have  used  the  following  Van-der-Waals-type  mixing  rules 


T^^  = (l-x)^T^^  + 2x(l-x)T^^^  + x^T^^, 


1(22) 


with  the  combining  rules 
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] 


3 


1(23) 


and 


1(24) 


where 


V = V Z ^/Z 
5 5 4 '5 


1(25) 


The  guides  in  establishing  the  parameter  values  of  k and  % are  the  experimental 
data  on  the  critical  line  and  vapor-liquid  equilibria.  Obtaining  the  location 
of  the  real  critical  line  from  the  surface  defined  by  1(2)  - 1(25)  is  not 


the  critical  line  is  approached  from  the  one-phase  region.  Since  the  critical 
temperatures  of  the  global  surfaces  for  isobutane  and  isopentane  are  about  2 K 
higher  than  the  physical  values,  we  required  for  the  critical  line  of  the  mixture 
surface  that  it  be  displaced  by  a similar  amount  from  the  measured  values.  The 
parameters  k and  2,  were  optimized  by  this  procedure  (Table  I 12) . 

The  global  surface  so  defined  for  the  mixture  was  used  in  the  comparisons 
with  the  critical-region  surface  developed  within  the  framework  of  this 
contract  (Part  II). 

2.3  A scaled  thermodynamic  surface  for  the  critical  region  of  isobutane 

The  global  thermodynamic  surface  that  we  developed  for  isobutane  is 
inaccurate  in  the  immediate  vicinity  of  the  critical  point.  In  a region  of 
the  size 


401.8  K < T < 414  K 


3 


173.4  P 322  kg/m 


and 


1(26) 


this  surface  should  not  be  used.  The  modern  theory  of  critical  phenomena  has 
provided  us  with  the  tools  to  describe  the  critical  region  of  fluids.  The 


15 


principle  of  critical-point  universality  implies  that  the  critical  anomalies 

of  all  fluids  and  fluid  mixtures  can  be  described  by  the  same  scaled  thermo- 

32  33 

dynamic  potential  save  for  two  arbitrary  scaling  constants.  ’ The 

34 

renormalization-group  approach  initiated  by  Kadanoff  and  formalized  by 
35 

Wilson  predicts  the  critical  exponents  and  amplitude  ratios  of  the  scaled 

36 

potential,  and  the  form  of  the  corrections  to  scaling  . This  potential  has 

37 

been  applied  to  a number  of  fluids,  so  that  a brief  summary  may  suffice 
here.  For  some  of  the  details,  we  refer  to  Part  II  of  this  report. 

The  scaled  thermodynamic  surface  is  a relationship  between  the  intensive 
thermodynamic  variables  pressure  P,  chemical  potential  p and  temperature  T. 

We  use  the  reduced  variables 

p=^. 


= H.  _ PP*  i 


T 


1(27) 


where  the  "starred"  and  "barred"  variables  are  those  defined  in  Sec.  2.1. 

The  fundamental  equation  yields  the  thermodynamic  potential  P as  a 
function  of  y and  T and  has  the  form 


P = Pq  (T)  + Ay  + P^^  Ay  AT  + AP 
with  AT  = T 4-  1 r 

Ay  y - PqCT)  . 


1(28) 


1(29) 


The  functions  Pq  (T)  and  y^  (T)  are  represented  by  truncated  Taylor  series 
expansions 

...  3 . . 

^0  (T)  = y^  4-  y.  (AT)^ 
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and 


1(30) 


P (T)  = 1 + Z P. (AT)^  . 

° j=l  ^ 

The  singular  part  AP  in  1(28)  is  related  to  Ay  and  AT  by  means  of  two 
auxiliary  (parametric)  variables  r and  0 

Ay  - r^‘^a0(l-e^)  , 

AT  - r(l-b^0^)  - cAy  , 


AP 


ak  p (0) 
0*^0 


^ ^3(6+l)+A^ 


akiPi(0) 


1(31) 


Here  8,  5 and  A^  are  universal  critical  exponents,  b a universal  parameter 

and  a,  c,  k and  k,  are  constants  that  differ  from  substance  to  substance, 
o 1 

2 

The  functions  Pq(9)  and  p2^(Q)  are  universal  quadratic  polynomials  in  0 given 
33  37 

elsewhere.  ’ The  parameters  of  the  revised  and  extended  scaled  potential 
P are  listed  in  Table  I 13.  Most  of  the  parameters  were  determined  by  a fit 
to  Beattie's  critical-region  PVT  data,  supplemented  by  data  points  generated 
from  the  global  surface  in  the  region  where  the  latter  is  still  valid,  but 
within  the  range  of  validity  of  the  scaled  equation.  The  expansion  coefficients 
\i2  and  y^  of  the  chemical  potential  were  obtained  by  fitting  to  c^  data  generated 
from  the  global  surface  in  the  same  range.  The  coefficients  y^  and  y^  were 
obtained  by  setting  the  energy  and  entropy  of  the  critical-region  surface 
equal  to  those  of  the  global  surface  at  the  reference  point  of  438  K in 

3 

temperature  and  205  kg/m  in  density. 

Detailed  information  on  this  surface  and  tables  of  thermodynamic  properties 
generated  from  it  can  be  found  in  Ref.  7. 

The  critical-region  surface  matches  with  the  global  surface  in  the 
supercritical  region  very  well.  There  are  some  problems  in  the  subcritical 
region  near  the  saturation  boundary.  There  are  very  few  reliable  data  in  that 
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region;  also  the  ranges  of  validity  of  the  global  and  the  critical  surface 
do  not  overlap  below  T . ' 'i. 
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Part  II 


A Scaled  Fundamental  Equation  for  Mixtures  of  Isobutane 
and  Isopentane  near  the  Gas-Llquld  Critical  Line 

1.  Introduction 

In  Part  I,  Sec.  2.3  of  this  report  we  discussed  how  the  thermodynamic 
behavior  of  a pure  fluid  near  a critical  point  Is  represented  by  a non- 
classical  scaled  potential.  In  this  part  of  the  report,  we  explore  the 
application  of  the  scaling  laws  to  fluid  mixtures.  In  principle,  this  Is 
straightforward:  critical-point  universality  Implies  that  the  properties  of 
all  systems  In  the  class  of  Islng-llke  system  scale  In  the  same  way.  This 
class  encompasses  not  only  fluids  and  fluid  mixtures,  but  also  many  other 
systems,  such  as  uniaxial  ferromagnets  passing  through  their  Curie  point  and 
binary  alloys  undergoing  an  order-disorder  phase  transition.  The  general- 
ization from  one-component  to  two-component  fluid  mixtures  was  made  by 

OQ 

Griffiths  and  Wheeler  , who  pointed  out  that  the  mixture,  at  constant 

OQ 

activity,  scales  like  the  pure  fluid.  Leung  and  Griffiths  subsequently 

developed  a model  based  on  this  principle,  and  applied  it  to  the  mixture  of 
3 . 4 . 

helium  and  helium.  For  the  underlying  model  for  the  pure  fluid,  they 
assumed  the  synanetry  of  the  simple  lattice  gas,  and  used  apparent, 
rather  than  asymptotic  values  of  the  critical  exponents,  as  an  approximate 
way  of  including  corrections  to  scaling.  Their  model  gave  an  acceptable 
description  of  near— critical  helium  mixtures  in  roughly  the  same  range 
of  reduced  variables  as  the  scaled  equations  we  discussed  in  Part  I. 

Holdover,  Rainwater  and  coworkers^®’  then  modified  the 

Leung-Grif f iths  model  in  order  to  obtain  more  accurate  results  for 
vapor-liquid  equilibria  in  hydrocarbon  mixtures  with  quite  dissimilar 
molecules.  They  introduced  terms  that  break  the  particle  - hole  symmetry 


and  used  aspects  of  the  law  of  corresponding  states  to  make  the  model  j 

more  predictive.  Thus  far,  their  model  has  not  been  applied  to  supercritical  j ; 

I 

mixtures.  I 

42  i ' 

Chang  and  Doiron  used  a minor  modification  of  the  original  Leung- 

Griffiths  model  to  describe  supercritical  carbon  dioxide  - ethane  mixtures. 

In  this  chapter  of  the  report,  we  first  present  a model  shaped 

after  that  of  Chang  and  Doiron  and  apply  it  to  isobutane  - isopentane  | 

mixtures.  Since  the  original  Leung-Grif f iths  model  cannot  handle  the 

refinements  of  liquid-vapor  asymmetry  and  corrections  to  scaling,  we 

first  had  to  construct  simple  scaled  potentials  (akin  to  that  of  a i 

lattice  gas)  for  pure  isobutane  and  isopentane  as  a starting  point  for 

constructing  the  thermodynamic  potential  of  the  mixture;  following  Leung 

and  Griffiths,  we  use  apparent  values  of  the  critical  exponents.  In 

Sec.  2 of  this  part  of  the  report,  we  develop  these  thermodynamic  potentials 

for  the  pure  fluids.  In  Sec.  3,  we  develop  the  Leung-Grif f iths  model  I 

and  in  Sec.  4 we  compare  the  model  with  experimental  data.  f 

After  this  report  was  completed,  we  discovered  a fundamental  difficulty 

I 

with  our  version  of  the  Leung-Grif f iths  model,  a difficulty  not  present 
in  the  original  model  or  the  versions  of  Holdover  et  al. , or  Chang  and 

Doiron.  This  difficulty  is  discussed  in  Section  4.3.  jj 

1 

In  the  latter  part  of  this  chapter,  we  present  a model  shaped  after  ’ 

I 

that  of  Holdover,  Rainwater  et  al. , and  show  that  it  represents  the  VLE  ; 

i 

properties  of  the  mixture  as  well  as  they  are  known  from  the  global  jl 

surface  and  the  available  data.  We  conclude  that  this  version  of  the  ! 

Leung-Grif f iths  model  is  capable  of  representing  the  two-phase  data  and 
that  there  is  therefore  no  reason  why  the  original  model  should  net  work 
reasonably  well  in  both  the  one  - and  two-phase  region. 

Since  this  part  of  the  report  was  developed  independently  from  the 
other  parts,  there  are  some  redundancies  and  some  differences  in  notation; 
eliminating  these  would  have  reduced  the  clarity  and  readability  of  this 
part  of  the  report.  


2.  Scaled  fundamental  equation  for  one-component  fluids 


2.1  Critical  power  laws 

Let  T be  the  temperature,  p the  density,  P the  pressure,  V the  volume,  y 
the  chemical  potential  per  mole,  A the  Helmholtz  free  energy,  S the  entropy,  H 
the  enthalpy,  the  heat  capacity  at  constant  volume,  Cp  the  heat  capacity  at 
constant  pressure  and  Kj  s p ^(3p/8P)^  the  isothermal  compressibility.  In 
practice  we  find  it  convenient  to  introduce  a symmetrized  compressibility 
defined  as 


X'j 


11(1) 


The  thermodynamic  properties  are  made  dimensionless  by  expressing  them  in 
terms  of  the  critical  temperature  T^,  the  critical  density  p^  and  the  critical 
pressure  P^,. 

T p P T 

X S.  0 p = 

T » P p ’ T P * 

c c 


y p T U 1 AT 

H = — • — £_£  u = _ . A • — ^ 

Mt  P’  VP’  VTP 
c c c 


11(2) 
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21 


1 


^ ^ ^ ^ 

Note  that  the  reduced  extensive  properties  A,  U,  H,  S,  are  all  taken 

per  unit  volume  rather  than  per  unit  mass.  The  reason  is  that  the  singular 
part  of  the  extensive  thermodynamic  properties  per  unit  volume  appears  to  be 
an  approximately  symmetric  or  antisymmetric  function  of  p ” shovm  in 

earlier  papers^"^ * . As  a consequence  x-j  is  also  a symmetric  function  of 


'2^ 

p - p^,  while  the  compressibility  X^/P  is  not. 

introduce  the  reduced  differences 

In  addition  we 

AT  = T + 1 , 

II(3a) 

Ap  = p - 1 , 

^ ^ 

II(3b) 

Ay  = y - y(p^,T)  . 

II(3c) 

To  represent  the  singular  thermodynamic  behavior  of  fluids  in  the  vicin- 
ity of  the  critical  point,  we  first  define  critical  power  laws.  The  exponents  i 
of  these  power  laws  depend  on  the  property  considered  and  on  the  path  along 
which  the  critical  point  is  approached.  The  most  important  critical  power  ' 

laws  are 


Cy  1aT1"“  (Ap  = 0.  AT  > 0)  , II(4a) 
Ap  = ± B|A^j^  (coexistence  boundary)  , II(4b) 
\ = rlATl"^  (aS  = 0,  at  > 0)  , II(4c) 
Ay  = ± D|Ap!^  (AT  = 0)  . II  (4d) 
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I The  critical  exponents  a,  3>  Y»  « are  related  by  the  exponent  relations 

I 

3(6+1)  = 2-a  , 3(6-1)  = y . 11(5) 

The  values  of  the  critical  exponents  are  universal,  i.e.,  they  are  the  same  for 

I 

all  systems  belonging  to  the  same  universality  class.  Theory  predicts  that 
fluids  near  the  gas-liquid  critical  point  belong  to  the  universality  class  of 
3-dimensional  Ising-like  systems  , The  exponents  for  this  universality  class 
are  currently  known  with  considerable  accuracy . It  is  noted  that  the 
critical  power  laws  defined  by  11(4)  are  the  leading  terms  of  expansions  around 
the  critical  point 

2,2  Linear-model  fundamental  equation 

As  our  fundamental  equation  we  consider  the  potential  P as  a function  of 
T and  p.  This  potential  satisfies  the  differential  equation 

dP  = UdT  + pdy  . 11(6) 

To  formulate  a scaled  fundamental  equation  for  the  thermodynamic  prop- 
erties of  fluids  near  the  critical  point  two  procedures  are  currently 
available.  The  oldest  approach  is  to  formulate  a simple  scaled  equation  that 
incorporates  the  critical  power  laws  defined  in  11(4)  without  nonanalytic 
corrections  to  the  power  laws.  We  have  found  that  this  approach  yields  a 
reasonably  accurate  representation  of  the  thermodynamic  properties  of  fluids 
in  the  critical  region,  but  with  effective  values^^  for  the  critical  exponents 
that  differ  slightly  from  the  theoretical  values^^*^^  calculated  for  the  uni- 
versality class  of  Ising-like  systems,  A more  rigorous  approach  is  to  formulate 
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a revised  and  extended  scaled  fundamental  equation  that  incorporates  the  most 
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important  correction  terms  to  the  asymptotic  critical  power  laws  * , This 

revised  and  extended  scaled  fundamental  equation  has  been  shown  to  yield  a 
good  representation  of  the  thermod5niamic  properties  of  various  fluids  in  the 
critical  region  including  isobutane  . The  Leung-Grif f iths  model  is  a 
generalization  of  the  simple  scaled  fundamental  equation  to  mixtures.  A gen- 
eralization of  the  revised  and  extended  scaled  fundamental  equation  to 
mixtures  is  not  yet  available.  Before  we  are  able  to  apply  the  Leung- 
Grif  f iths  model  to  isobutane-isopentane  mixtures,  we  first  need  to  represent 
the  thermodynamic  properties  of  the  individual  components  near  their  critical 
points  in  terms  of  a simple  scaled  fundamental  equation. 

The  potential  P is  written  as  the  sum  of  a regular  part  P^  and  a scaled 
part  P^ 


P = P^(AT,Am)  + Pg(AT,Ap)  . 11(7) 


The  functions  P^(AT,Ay)  in  11(7)  and  y(p^,T)  = yg(AT)  in  II(3c)  are  analytic 
functions  that  can  be  approximated  by  truncated  Taylor-series  expansions: 


3 

y(p  ,T)  =I^(AT)  = y + s y.  (AT)^  , 11(8) 

c ^ c . , 1 


3 

P (AT, Ay)  = Pq(AT)  + Ay  = 1 + Z P.  (AT)^  + Ay  . 11(9) 

i=l 
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In  order  to  specify  the  scaled  part  of  the  potential  one  introduces  a new 
set  of  variables  r and  9 in  the  one-phase  region  that  are  related  to  the 
physical  variables  AT  and  Ay  by  the  transformation 


AT  = r(l-b^e^)  , Il(lOa) 

Ay  = ar^'^0(l-0^)  , Il(lOb) 

where  a and  b are  constants.  This  transformation  is  chosen  such  that  0=0 
corresponds  to  the  critical  isochore  Ap  = 0 above  the  critical  temperature, 

9 = ± 1/b  to  the  critical  isotherm  AT  = 0 and  9 = ± 1 to  the  two  branches  of 
the  coexistence  curve.  In  the  linear-model  scaled  equation  adopted  here,  the 

constant  b in  the  transformation  equation  II (10a)  is  related  to  the  critical 
, 44 

exponents  by 


2 _ 6- 

(6-l)(l-26) 


11(11) 


In  terms  of  the  parametric  variables  the  scaled  part  of  the  potential  is 
given  by 


P 

s 


^^6(6  + 1)  J g 


11(12) 


with 
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II(13a) 


) = 3<S-33-b  OY 

° 2b^(2-a)(l-a)a 


y 


p 


2 


e6-3g-b^a(2B6-l) 

2b^(l-a)a 


II(13b) 


236-3 

2a 


II(13c) 


With  the  potential  P specified  by  11(7),  11(9)  and  11(12)  all  thermodynamics 
properties  can  be  derived  with  the  aid  of  the  differential  relation  11(6). 
Specifically  one  obtains  for  the  density  the  simple  relation 


Ap 


= kr^e 


11(14) 


The  simple  scaled  fundamental  equation  adopted  here  is  known  as  the  "linear 
model"  because  of  the  linear  relationship  between  Ap  and  9.  A complete  list 
of  the  relevant  linear  model  equations  is  presented  in  Appendix  B, 

The  amplitudes  A,  B,  T,  D of  the  critical  power  laws  defined  in  11(4)  are 
related  to  the  linear-model  fundamental  equation  by"^^ 


A = ak(2-a)(l-a)p  , II(15a) 

o 

B = k/(b"-l)^  , II(15b) 

r = k/a  , II(15c) 

D = a(b^-l)b'^"^/k'^  . II(15d) 
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In  summary,  the  linear-model  fundamental  equation  contains  the  following 
parameters:  the  two  critical  exponents  3 and  6,  the  two  parameters  a and  k 

which  together  with  the  exponents  0 and  6 specify  the  amplitudes  of  the  crit- 
ical power  laws,  the  critical  parameters  T^,  p^,  P^,  three  pressure  background 
parameters  , P2 » four  caloric  background  parameters  y^,  y^  y2 

^ /s.  rv 

and  y^*  The  parameters  y^  and  y^^  specify  the  reference  points  of  energy  and 
entropy.  In  this  report  we  adopt  the  convention  that  the  energy  U and  the 
entropy  S of  isobutane  and  of  isopentane  are  taken  to  be  zero  at  their  respec- 
tive critical  points.  This  condition  implies 


y = 1 , II(16a) 

c 

y^  = P^  . II(16b) 

In  our  previous  work  we  have  found  that  the  simple  linear-model  funda- 
mental equation  can  be  used  in  a range  around  the  critical  point  approximately 
specified  by"^^. 


0.0005  < [at]  < 0.03  , |Ap|  < 0.25  , 11(17) 


if  for  the  critical  exponents  the  following  effective  universal  values  are  used 


B = 0.355  , 6 = 4.352  , a = 0.100  , y = 1.190  11(18) 

and  with  b as  determined  from  11(11) 

b^  = 1.3909  . 11(19) 
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2,3  Determinaton  of  parameters  for  isobutane 


The  scaled  fundamental  equation  presented  here  for  isobutane  is  valid  in 
the  range 


407.6  K < T < 420.1  K 


11(20) 


169  kg/m^  < p < 282  kg/m^ 


The  critical  temperature  and  the  critical  density  of  isobutane  have  been 
determined  experimentally  as^’^’^(c.f.  Table  I 11) 


T = 407.84  K 
c 


y 


11(21) 


= 225.5  kg/m^  = 3.8796  mol/dm^ 


In  this  part  of  the  report  we  adopt  the  convention  that  quantities  with  a 
superscript  (4)  refer  to  pure  isobutane  and  those  with  a superscript  (5)  to 
pure  isopentane. 

The  critical  pressure  and  the  parameters  a,  k,  , P^ , and  P^  can  be 
determined  by  fitting  the  linear-model  equation  to  experimental  equation-of- 
state  data.  The  primary  experimental  data  available  for  this  purpose  are  the 
data  reported  by  Beattie  et  al.^^’^^  Unfortunately  the  experimental  data  of 
Beattie  et  al.  cover  only  a small  portion  of  the  range  specified  by  11(20). 
Specifically,  we  have  174  data  points  restricted  to  temperatures  between 
407.764  K and  408.314  K and  densities  between  180,5  kg/m^  and  270.2  kg/m^,To 
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cover  the  entire  range  specified  in  11(20),  we  supplemented  the  experimental 

data  with  pressure  data  calculated  from  the  revised  and  extended  scaled 

equation  earlier  constructed  for  isobutane  by  Levelt  Sengers  et  al.^ 

Specifically,  we  calculated  24  data  points  spaced  by  5 kg/m  intervals  in  the 

3 3 

density  range  168  kg/m  < p < 283  kg/m  at  each  of  the  following  temperatures: 
412  K,  414  K,  416  K,  418  K and  420  K.  We  thus  obtained  120  calculated  data 
points  supplementing  the  174  experimental  data  of  Beattie  et  al.  We  assigned 
to  the  experimental  and  to  the  calculated  data  points  the  error  estimates  of 
Beattie  et  al.,  namely 

Op  = 10"“^  MPa  , = 0.005  K , = 0.058  kg/m^  . 11(22) 

It  was  found  that  the  linear-model  equation  represents  these  data  with  a 
reduced  chi-square  of  1.84.  The  parameter  values  deduced  from  the  fit  are 
given  in  Table  II  1. 

The  parameters  and  y^  were  selected  in  accordance  with  11(16).  The 
parameters  ^2  and  y^  can  be  determined  from  experimental  specific  heat  data. 

In  the  absence  of  experimental  specific  heat  data  for  isobutane  in  the  crit- 
ical region,  we  resorted  again  to  the  revised  and  extended  scaled  equation  for 
isobutane  from  which  we  generated  calculated  data  for  the  specific  heat  at 
constant  volume.  A total  of  168  points,  at  seven  temperatures  at  2 K intervals 
between  408  K and  420  K and  at  5 kg/m  density  intervals  in  the  range  from 
168  kg/m"^  to  283  kg/m^,  were  used  in  the  fit  to  determine  ^2  the  values 

thus  obtained  for  these  parameters  are  given  in  Table  II  1.  The  differences 
between  values  calculated  from  the  simple  scaled  equation  presented  here  and 
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those  calculated  from  the  original  revised  and  extended  scaled  equation  are 
about  3%  very  close  to  the  critical  point  and  are  well  within  1%  at  the  edges  of 
the  range  specificed  by  11(20).  The  differences  between  Cp  values  calculated 
from  the  simple  scaled  equation  and  those  calculated  from  the  revised  and 
extended  scaled  equation,  are  generally  about  3%  in  the  range  11(20),  but  very 
close  to  the  critical  point,  specifically  on  the  isotherm  at  408  K,  the 
differences  become  as  large  as  20%. 

2.4  Determination  of  parameters  for  isopentane 

The  scaled  fundamental  equation  presented  here  for  isopentane  is  valid  in 
the  range 


460.3  K < T < 474.3  K 


177  kg/m^  < p ^ 294  kg/m"^ 


11(23) 


The  critical-point  parameters  of  isopentane  are  (c.f.  Table  I 11): 


T (5)  = 
c 


P 

c 


(5)  = 


460.51  K 


235.57  kg/m^  = 


3.371  MPa 


3.265  raol/dm' 


11(24) 
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In  the  absence  of  experimental  PVT  data  for  Isopentane  in  the  critical 
region,  we  cannot  determine  the  parameters  a and  k in  the  scaling  functions 
from  experimental  data  directly  although  the  product  ak  can  be  inferred  from 
vapor  pressure  data.  Instead  we  assume  that  the  singular  parts  of  the  poten- 
tials of  isobutane  and  isopentane  satisfy  corresponding  states  which  implies 
a and  k to  be  the  same  for  isobutane  and  isopentane.  Previous  work  has  shown 
this  assumption  to  be  valid  when  the  thermodynamic  properties  of  D2O  in  the 
critical  regiona  are  compared  with  those  of  H2O  in  the  critical  region,  even 
though  these  fluids  themselves  do  not  satisfy  corresponding  states^*^. 

We  determine  the  pressure  background  parameters  Pj^,  ^3  from  exper- 

imental vapor-pressure  data  given  in  Table  I 6.  The  vapor-pressure  data  are 
represented  by  the  simple  scaled  equation 

P = 1+  P^AT  + P^  + P2(AT)^  + P^CAT)^  , 11(25) 

where  the  coefficient  P^  of  the  nonanalytic  term  is  (c.f.(3.20)  in  Appendix  B) 


akCp^+p^+p^) 
’’a  ' (j_.„2^6(o+l) 


11(26) 


With  the  previously  adopted  values  of  b^,  a and  k we  find  P^  = 30.2137.  With 
estimated  experimental  errors 


= 0.01  K , Op  = 0.003  MPa 


11(27) 
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Che  vapor-pressure  equation  11(25)  reproduces  Che  experimental  data  with  a i 

chi-square  of  0.52  and  with  coefficients 

I 

= 6.1720,  P2^^^  = -18.0281,  ^3^^^  = 19.4960  . 11(28)  j 

A detailed  comparison  between  the  experimental  and  calculated  vapor  pressures  | 

i 

is  presented  in  Table  I 6.  | 

/V  [ 

In  order  to  determine  the  caloric  background  parameters  \i2  ^3  need 

estimates  for  the  specific  heat.  Again  experimental  information  for  the 

I 

specific  heat  of  isopentane  in  the  critical  region  is  not  available.  As  an  i 

alternative  we  decompose  the  specific  beat  into  the  ideal-gas  specific  heat  j 

i 

C,,  .j  and  an  excess  specific  heat  AC  = C - C . , | 

V ,id  ^ V V V ,id  1 


C = AC  + C , , , 

V V v,id 


11(29) 


and  assume  that  the  excess  specific  heat  contributions  AC  and  AC  for 

v V 


isobutane  and  isopentane  satisfy  corresponding  states.  For  isobutane  the  ideal- 


,(4) 


gas  specific  heat  C^  is  calculated  from  the  global  thermodynamic  surface  of 


3 (4)  I 

Waxman  and  Gallagher  . The  excess  specific  heat  AC  for  isobutane  is  then  ;! 


deduced  as  C , with  again  calculated  from  the  revised  and 

V V ,id  V “ 


extended  scaled  equation  of  state  for  isobutane  given  in  Part  I.  The  excess 


specific  heat  [ “(5) » ~(5) ) isopentane  is  then  constructed  from 


.(4)/t 


4C''  7 — , •^//\\  for  isobutane  as 


V (4) 


(4)' 


AC  (5) 


AC 

V 


11(30) 
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with 


Z 

c 


(4) 


p (4) 

R p 


p (5) 


c c 


11(31) 


Finally,  the  total  specific  heat  of  isopentane  is  obtained  as  the 

sum  of  the  excess  specific  heat  estimated  on  the  basis  of  corresponding 

states  and  the  ideal-gas  specific  heat  of  isopentane  calculated  from  the 

generalized  corresponding  states  surface  described  in  Part  I,  Sec.  2,2.  Having 
thus  generated  a set  of  data  for  isopentane  we  fit  the  simple  scaled  equaton 
to  these  data  to  determine  ^2  ^3*  3’he  differences  between  the  generated 

specific  heat  data  and  the  values  calculated  from  the  simple  scaled  equation  are 
at  most  2.2%  very  near  the  critical  point  and  within  1%  in  most  of  the  range 


specified  by  11(23),  The  coefficients  and  again  are  fixed  by  the 
convention  that  the  energy  and  entropy  of  isopentane  are  taken  to  be  zero  at 
the  critical  point. 

The  parameters  for  isopentane  are  all  given  in  Table  II  1. 


3.  Scaled  fundamental  equation  for  fluid  mixtures 


3.1  Reduced  thermodynamic  variables 


The  scaled  fundamental  equation  for  fluids  as  discussed  in  Part  II,  Sec.  2 

39 

can  be  generalized  by  a procedure  proposed  by  Leung  and  Griffiths  . However, 
the  equations  used  by  Leung  and  Griffiths'^  as  well  as  by  others^  are  not 
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in  dimensionless  form  unlike  the  scaled  equations  used  in  this  report. 


This 


makes  the  equations  inconvenient  when  one  wants  to  change  from  one  system  of 
units  to  another.  As  a first  step,  therefore,  we  reformulate  the  equations  of 
Leung  and  Griffiths  so  as  to  bring  them  into  dimensionless  form. 

We  treat  pure  isobutane  as  the  reference  fluid  and  we  make  the  thermo- 
dynamic properties  of  the  mixtures  dimensionless  with  the  aid  of  the  critical- 

(4)  (4)  (4) 

point  parameters  and  of  isobutane,  which  are  the  "starred” 

reference  values  listed  in  Table  I 9.  We  thus  define  for  the  mixture 


T 


U 


U 


V 


11(32) 


Let  and  p^  be  the  m.olar  densities  and  and  the  chemical  potentials  of 
isobutane  and  isopentane  in  the  mixture.  These  quantities  are  similarly  made 
dimensionless  as 


11(33) 


4 c c 
7"  . (4) 


P,  P ^^^T 
5 c c 

7"  "TTT 


11(34) 
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'if  In  this  report  we  also  consider  the  mole  fraction  x of  isopentane  in  the 


mixtures  which  is  related  to  the  densities  and  by 


X = 


P4  + P5 


11(35) 


As  in  the  case  of  the  one-component  fluids  we  shall  formulate  a funda- 
mental  equation  for  the  potential  P.  This  potential  P satisfies  the 
differential  equation 


dP 


p5dlJ5  UdT 


11(36) 


OQ  /V 

Following  Leung  and  Griffiths  we  do  not  consider  P as  a function  of 
and  T,  but  as  a function  of  a new  set  of  variables  4,  x and  h which  are 
related  to  the  original  set  of  intensive  variables  by 


P4 

e + e ‘ 


11(37) 


T 


T - T^(4) 


11(38) 


^4  P5 

ln(e  + e ) - H(c,t) 


11(39) 
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The  variable  ^ plays  the  role  of  an  intensive  composition  variable  such  that 
5=0  corresponds  to  pure  isobutane  and  5 = 1 to  pure  isopentane.  The  critical 

t 

temperature  T^(5),  the  critical  density  P^(C)  and  the  critical  pressure 
are  now  treated  as  functions  of  5 rather  than  of  the  mole  fraction  x.  The 
critical  parameters  of  the  mixture  in  dimensionless  form  are 


T (4) 

T (?)  ^ 

T^(?) 


Pc(?) 

=— TO 


P (O 

T (5) 


11(40) 


^4 

The  function  H(5,t)  represents  the  value  of  ln(e  + e ) along  the  phase 
boundary.  The  variables  x and  h are  defined  in  such  a way  that  in  the  limit 
5 ->■  0 they  reduce  to  the  variables  AT  and  Ay  for  pure  isobutane. 

3.2  Scaled  fundamental  equation  for  the  mixture 

In  analogy  to  11(7)  we  decompose  the  potential  P again  in  a regular  part 
P^  and  a scaled  part  P^ 


P = P^(5,T,h)  + P^(5,T,h) 


11(41) 


The  functions  P^(5,T,h)  in  11(41)  and  H(5,t)  in  11(39)  are  analytic  functions 
that  are  represented  by  truncated  Taylor  series  expansions  in  analogy  to  11(8) 
and  11(9): 


3 

H(5,t)  = Z H (5)t^  , 11(42) 

i=0 
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11(43) 


P,(C.T,h)  = Z P.(Ot^  + P^(Oh 
1=0 


The  singular  part  P^  is  specified  as 

Pg(4,T,h)  = q(5)AP(T,h) 


11(44) 


with 


T = i(5)T  . 11(45) 

The  function  AP(T,h)  is  to  be  identified  with  the  function  AP(A’?,Ap)  defined 
in  Appendix  B for  one-component  fluids.  Specifically 

T = Jl(c)T  = r(l-b^0^)  , 11(46) 

h = ar^'^Sd-e^)  , 11(47) 

AP  = akr®^^'^^^(pQ  + p^e^  + p^8^)  . 11(48) 

The  functions  T^(c),  p^(?),  P^(?)  = Pq(4),  \(C),  ^^(4).  q(C).  Ji(C)  are  all 
analytic  functions  of  the  variable  C and  may  be  expanded  in  series  of  the  form 

f(?)  = fQ+f^C+f^dl-O  + f34^(l-4)  + ...  . 11(49) 
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We  note  that  the  coefficients  a and  k in  11(47)  and  11(48)  are,  in  principle, 
functions  of  the  variable  z,.  However,  as  discussed  in  Part  II,  Sec.  2.4,  for 
our  mixture  a = a^^^  = a^^^  and  k = k^^^  ==  k^^^  are  taken  as  constants, 

A list  of  equations  for  various  thermodynamic  properties  implied  by  this 
fundamental  equation  is  presented  in  Appendix  C. 

3,3  Determination  of  parameters  for  isobutane-isopentane  mixtures 

The  equations  presented  in  Part  II,  Sec.  3.2,  specify  the  thermodynamic 
surface  of  the  mixture  in  terras  of  the  following  constants  (j  = 0,  1,  2,  3): 


In  determining  these  constants  we  need  to  impose  the  conditions  Chat  the 
fundamental  equation  reproduce  the  fundamental  equation  for  pure  isobutane  in 
Che  limit  C 0 and  the  one  for  pure  isopentane  in  the  limit  C ^ This 
condition  implies  the  following  relationships  between  Che  constants  in  the 
fundamental  equation  of  the  mixture  and  chose  in  Che  fundamental  equation  of 
the  pure  components: 


T 

c 

T 

c 


11(50) 


^cl 


D 


1 


11(51) 


1 


11(52) 


p (5)  ^ (4) 
^ c c 

^00  " ^ ’ ^01  " TTsT 


p 


iO 


^00 


7 (5) 

, „ c - (5) 

bl  = 7W 


> 


q 


0 


= 1 , 


i 


0 


= 1 , 


11(53) 


11(54) 


11(55) 


11(56) 


11(57) 
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The  first  two  terms  in  the  expansion  11(49)  for  T^C?),  p^(^),  P^C^), 

9(4) > ^(?)  are  thus  completely  specified  by  the  properties  of  the  pure 

components.  In  order  to  determine  the  coefficients  of  any  higher-order  terms 

in  the  expansions  we  need  thermodynamic  property  data  for  the  mixture.  For 

this  purpose  we  consider  critical  line-data  and  data  for  the  dew-bubble  curve 

discussed  in  Part  I and  listed  in  Tables  I 5 and  I 7. 

The  available  experimental  information  is  not  sufficiently  extensive  to 

determine  coefficients  of  the  higher  order  terms  in  the  expansion  11(49)  for 

all  of  the  functions  T (5),  p (^),  P, (5),  H (5),  q(c),  £(?).  In  practice  we 

cell 

only  retain  higher  order  terms  in  the  expansions  for  p^(c),  Pq(?)  and  Pj^(c) 
and  we  restrict  ourselves  to  the  following  truncated  expansions: 


11(58) 


11(59) 


11(60) 


11(61) 


11(62) 


11(63) 
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Hi(0  = , 

q(?)  = ^0 

1(0  = ^0 


11(64) 


11(65) 


11(66) 


We  repeat  that  the  coefficients  of  the  constant  and  linear  terms  in  these 
expansions  are  determined  by  the  constants  in  the  equation  for  pure  isobutane 
and  isopentane.  The  coefficients  of  the  higher  order  terms,  i.e.  P^2*  ^02’ 
Pq2  ^nd  Pj^2  3re  estimated,  primarily  on  the  basis  of  the  critical-line  data 
presented  in  Table  I 5,  as 


p^2  " 0.06 


^02  » ^03  ""  -0-03 


11(67) 


P 


12 


= 7 


A complete  set  of  coefficients  for  isobutane-isopentane  mixtures  is 
presented  in  Appendix  D.  Tables  of  thermodynamic  properties  calculated  from 
the  scaled  fundamental  equation  are  presented  in  Appendix  E, 
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4,  Discussion  of  results 


4,1  Comparison  with  experimental  data 

The  scaled  fundamental  equation  for  the  mixture  was  constructed  in  such  a 
way  as  to  reproduce  the  thermodynamic  properties  of  the  pure  components.  In 
addition  we  can  make  a comparison  with  the  experimental  data  for  the  critical 
point  and  the  coexistence  curve  data  for  the  mixture  as  a function  of  concen- 
tration presented  in  Tables  I 5 and  17. 

The  critical  parameters  of  the  mixtures  are  given  by  the  functions  T^C^), 
p^C?)  and  where  is  related  to  the  mole  fraction  x of  isopentane  by 


p(0 


11(68) 


On  the  critical  line  p(c)  = P (?)  and  Q(?)  = Q (?)  with 

c c 


Q,(0  = 


d? 


- P^(?> 


dT  (O 
c 

d? 


- P (?) 
c 


"dH  (O  dT  (O 


11(69) 


In  Taole  II  2 we  list  the  critical  parameters  of  the  mixture  calculated  from 
this  equation  for  various  values  of  the  composition;  these  calculated  values 
are  to  be  compared  with  the  experimental  values  given  earlier  in  Table  I 4, 
The  comparison  is  made  in  Fig,  II  1.  We  consider  the  agreement  with  the 
experimental  data  satisfactory. 
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The  experimental  coexistence  curve  data  presented  in  Table  I 7 are  compared 
with  coexistence  curve  data  calculated  from  our  fundamental  equation  at  constant 
compositions  x in  Fig.  II  2.  Here  the  agreement  is  only  approximate;  there  are 
systematic  disagreements  between  the  experimental  and  calculated  slopes  of  the 
coexistence  curve  data. 

We  tried  to  improve  the  agreement  between  experiment  and  theory  by  allowing 
for  nonzero  values  of  the  coefficients  Hq2  and  Hj^2  expansion  11(49) 

for  Hq(5)  and  This  procedure,  however,  did  not  lead  to  a significant 

change  in  the  values  obtained  for  the  thermodynamic  properties  from  our  funda- 
mental equation. 

4.2  Comparison  with  global  thermodynamic  surface  for  isobutane-isopentane 
mixtures 

The  application  of  the  scaled  fundamental  equation  to  the  fluid  mixture 
is  at  any  given  composition  restricted  to  a region  around  the  critical  point 
approximately  given  by  11(17).  In  Part  I we  have  developed  a global  analytical 
fundamental  equation  using  generalized  corresponding  states.  This  formulation  is 
valid  far  away  from  the  critical  line  but  ceases  to  be  valid  near  the  critical 
line.  One  can  investigate  the  consistency  of  the  two  procedures  by  comparing  the 
two  thermodynamic  surfaces  at  temperatures  and  densities  where  both  surfaces  are 
expected  to  be  applicable.  In  Part  I,  we  discussed  such  a test  for  pure 
isobutane  and  found  good  consistency  of  the  global  and  scaled  surfaces  at 
supercritical  conditions.  Here,  we  discuss  such  a comparison  for  an  equimolar 
mixture  of  isobutane-isopentane  (x  = 0.5);  at  this  composition  the  differences 
between  the  two  prediction  methods  are  likely  to  be  largest . 
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An  equimolar  mixture  of  isobutane  and  isopentane  has  a critical  temper- 
ature  of  437.71  K and  a critical  density  of  3.569  mol/dm'^.  From  11(17)  it 
follows  that  our  scaled  fundamental  equation  should  be  valid  at  temperatures 
up  to  452  K and  at  densities  from  2.7  mol/dm  to  4.5  mol/dm  . In  Table  II  3 
we  present  a comparison  between  the  pressure  and  its  first  derivatives  as 
calculated  from  the  global  analytical  fundamental  equation  at  452  K and  those 
calculated  from  the  scaled  fundamental  equation.  It  is  seen  that  the  differ- 
ences are  still  substantial.  The  differences  in  the  pressures  are  of  the  order 
of  5%  and  for  the  derivatives  they  increase  up  to  15%  for  the  thermal  pressure 
coefficient  and  up  to  20%  for  the  compressibility.  In  the  absence  of  direct 
experimental  PVT  data  for  this  mixture  one  cannot  ascertain  which  of  the  two 
surfaces  yields  the  more  accurate  values  at  this  temperature.  Moreover,  the 
procedure  for  applying  these  fundamental  equations  to  the  mixture  is  not 
unique.  If  we  compare  the  pressures  calculated  at  452  K from  our  fundamental 
equations  with  those  calculated  from  the  global  surface  but  with  a different 
set  of  coefficients,  we  find  differences  of  the  order  of  3%  but  with  the 
opposite  sign!  We  conclude  that  the  differences  quoted  in  Table  II  3 are  a 
measure  of  the  uncertainty  with  which  the  properties  of  these  mixtures  are 
predicted  with  currently  available  techniques,  when  the  available  data  base 
for  the  mixture  is  modest  to  scant. 

4.3  Concluding  remarks 

39 

The  procedure  devised  by  Leung  and  Griffiths  leads  to  a scaled  funda- 
mental equation  that  correctly  reproduces  the  scaled  fundamental  equation  for 
the  pure  components  near  their  respective  critical  points.  In  order  to 
calculate  the  properties  of  the  mixture  one  needs  to  specify  the  following 
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functions:  p^C^),  P^(?),  H^(c),  q(^)  and  Jl(c)  as  a function  of  For 

this  purpose  we  considered  an  expansion  of  the  form  (of  11(49)) 

f(c)  = fg  + + f24(l-C)  + f3C^(l-0  + ...  . 11(70) 

The  coefficients  fg  and  fj^  for  all  the  functions  mentioned  above  are  completely 

determined  by  the  properties  of  the  pure  components  as  shown  in  Sec.  3.3 
of  Part  II.  Therefore,  if  we  were  to  assume  that  all  functions  were  linear 
functions  of  ^ , we  could  predict  the  properties  of  the  mixtures  from  the 
properties  of  the  pure  components.  Unfortunately,  the  assumption  that  T^(5), 
P^(?)>  are  linear  functions  of  5 is  not  justified. 

In  our  approach  we  took  the  energy  and  the  entropy  of  the  pure  components 
to  be  zero  at  their  respective  critical  points.  This  choice  was  made  so  as  to 

bring  the  equations  into  their  most  simple  form.  A different  choice  of  the 

reference  points  for  zero  energy  and  entropy  would  affect  the  variable  C 
through  its  definition  11(37).  Hence  if  the  functions  f(c)  were  linear  for 
one  choice  of  reference  points  for  zero  energy  and  entropy,  they  would  no 
longer  be  linear  for  another  choice.  We  do  not  know  a priori  for  which  choice 
the  linear  approximation  would  be  optimum.  This  problem  has  also  been  men- 
tioned by  Leung  and  Griffiths  in  their  original  paper.  In  our  approach  we 
took  the  reference  points  as  the  critical  points  of  the  pure  fluids  and 
retained  nonlinear  terms  in  the  expansions  for  p^(c),  P^(^)  only. 

However,  if  we  a take  another  set  of  reference  points  for  zero  energy  and 
entropy  and  repeat  the  same  procedure  described  in  this  report  we  do  obtain 
different  values  for  the  pressures  of  the  mixture. 
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The  tables  presented  in  this  report  are  the  best  we  can  obtain  with  the 
information  and  methods  presently  available.  However,  in  order  to  develop  a 
method  which  can  truly  predict  the  thermodynamic  properties  of  fluid  mixtures 
near  the  critical  points,  it  will  be  necessary  to  develop  a procedure  for 
applying  the  scaling  laws  to  mixtures  in  a manner  which  is  less  sensitive  to 
the  choice  of  the  reference  points  of  zero  energy  and  entropy. 
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4.4  Critique 

I After  completion  of  the  work  reported  in  4.1  - 4.3,  a problem  was 

I discovered  with  the  definition  of  the  variable  i,  in  terms  of  the  reduced 

'1  Oi  'V 

I chemical  potentials  y^,  11(37).  At  x near  0,  by  virtue  of  thermodynamics 
« we  have  y^  ^ RT£n  x,  while  y^  RT  £n(l-x)  for  x near  1.  Since 
; ^ ^4-1 

ti  y^  = ^c4’  ^c4  critical  compressibility  factor  of  pure 

^ -P  -1 

I isopentane,  C will  be  proportional,  for  small  x,  to  x , with  p=Z^^, 

li  whereas  in  the  original  Leung-Grif f iths  model  and  in  the  variants  of 

Chang  and  Doiron  and  of  Holdover  et  al.  ? is  proportional  to  x for  small 

< X.  Although  this  nonanalytic  dependence  of  ? on  x is  not  forbidden  in 

i 

l!  principle,  it  is  inconsistent  with  the  asstimption  that  the  critical  line 
i parameters  and  the  background  functions  are  analytic  in  Thus,  the 

I 

way  the  Leung-Grif f iths  model  was  made  dimensionless  has  generated  a 

! 

' deeply  - hidden  inconsistency  that  is  very  likely  the  cause  of  the 

I 

I difficulties  we  experienced  with  the  model  and  that  will  have  to  be 

i 

removed. 
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5.  An  alternate  model  for  VLE  in  the  critical  region 

39 

We  now  present  an  alternate  version  of  the  Leung-Griffiths  model 

developed  by  Holdover,  Rainwater  and  co-workers^^’  which 

yields  a much  improved  representation  of  VLE  for  the  isobutane-isopentane 

mixture  in  the  critical  region.  Among  other  modifications,  this  alternate 

model  allows  for  a finite  slope  in  the  rectilinear  diameter  of  the 

temperature-density  coexistence  curve.  Consequently,  it  distorts  the 

line  of  symmetry  in  the  supercritical  one-phase  region,  and  therefore  is 

likely  to  present  difficulties  in  the  fitting  of  one-phase  thermodynamic 

mixture  data,  although  it  has  yet  to  be  tested  for  this  purpose.  On  the 

other  hand  the  model  has  provided  successful  VLE  fits  in  pressure, 

temperature,  and  density  for  a range  of  binary  mixtures-,  including  some  ; 
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with  relatively  wide  dew-bubble  curves  * which  are  hard,  tp_ fit  by 
other  methods. 

The  field  variable  space  u (5,  t,  h)  as  originally  defined  by  Leung 
and  Griffiths  is  retained,  where  R is  the  gas  constant  and 


RT  * 


u./RT 
e 5 


y,/RT  y./RT  * 
K e + e ^ 


11(71)  I 


11(72) 


T » (RT^(^))"-  - (RT)”^, 


11(73) 


_ y,,/RT  _ _ 

h = £n  (K  e + e ) -H  (?,  x)  , 


11(74) 


- - Pc/RT 

H (4,  t)  = iln  (K  e + e ^ )q. 


11(75) 


The  variables  x and  h have  been  written  with  bars  in  order  to 
distinguish  them  from  x and  h defined  by  Eqs.  11(37)  - 11(39).  The 


thermodynamic  potential  oj  is  proportional  to  P,  Eq.  11(32),  and  x is 
proportional  to  x,  Eq.  11(38),  whereas  the  exponents  in  the  definitions  of 
C and  h are  proportional  to  the  corresponding  exponents  in  the  definitions 
of  Eq.  11(37),  and  h,  Eq.  11(39).  However,  ^ and  h are  not  themselves 
proportional  to  ^ and  h.  K is  a constant,  in  effect  determined  by  the  choice 
of  zeroes  of  and  y^,  and  the  subscript  0 denotes  the  value  on  the 
coexistence  surface. 

As  the  mole  fraction  of  isopentane  goes  to  zero,  y^  diverges  as 

y^  'v  RT  £n  x 11(76) 

and  thus  for  dilute  solutions  of  isopentane  C is  linear  in  x,  as  it  should  be. 

The  thermodynamic  potential  has  been  explicitly  defined  in  Appendix 
A of  Ref.  51,  and  essentially  follows  Eqs.  11(42)  -11(48)  in  terms  of  the 
above  alternate  variables,  except  that  no  analog  of  il(C)  , Eq.  11(46),  exists 
and  that  x is  replaced  by  a new  variable  t, 

T - T^(0 

t = ^ , 11(77) 

lyo 

which,  in  contrast  to  X of  Eq.  11(38),  is  S3m3metric  on  interchange  of  the 
two  pure  fluid  labels. 

The  near-critical  coexistence  region  is  defined  by  0 _<  C ^ 1,  -0.1  < t < 0, 
and  h = 0.  In  the  present  method,  it  is  asstimed  that  the  vapor-pressure 
curves  and  coexisting  densities  of  the  pure  fluids,  as  well  as  the  critical 
line  (i.e.  T^  (x) , (x)  and  (x))  are  known.  The  objective  is  to  fit 

the  model  to  the  experimental  dew-bubble  curves  and  mixture  temperature-density 
coexistence  curves,  or,  equivalently,  to  specify  the  liquid  density  p^, 
vapor  density  p^,  liquid  composition  x,^  and  vapor  composition  x^  as 
functions  of  ^ and  t.  Unlike  the  model  of  Secs.  II.  3 and  II.  4,  chemical 
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I 

I 

I 

I 

1 

potentials  from  the  pure-fluid  equations  of  state  are  not  used  as  input.  j 

The  loci  of  constant  X are  given  by  the  equation  ' 

i 

T = 1 + 0,(5)  + c,(?)  t + tAX)  + CAi>  . II(  j 

T p (-J  3 4 5 6 

^ I 

where  ! 

(?)  = + ? (C<®)  - II  a 

for  i = 3,  4,  5s  6*  and  the  superscripts  refer  to  pure  isobutane  and  pure  ' 

(A)  fs)  I 

isopentane.  According  to  the  principle  of  corresponding  states,  • * 

The  above  equations  imply  that  lines  of  constant  ^ form  a^setrof-nearly,  paralld 

I 

curves  (cf.  Fig.  1 of  Ref.  41)  and  that  the  principle  of  corresponding  states  j 
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is  obeyed  approximately  by  the  mixture  in 'field  variable  space  . i 

The  values  of  and  are  listed  in  Table  II  4 where,  for  consistencj 
with  Eqs.  Il(lOa)  and  11(12),  i 

I 

= a k (Pq  + pg  + p^)!!  - = 30.22.  ^ 11(0 

The  remaining  coefficients  are  found  from  a linear  least  squares  fit. of  ^ : - | 

(P/T).  {IJV)  - Co  (~t)^"°^  -1  to  a cubic  polynomial  by  means  of  the  data  of  ' 
c c j I 

Tables  I 3 and  I 6. 

i 

The  coexisting  densities  are  given  by  ' r :r  , :| 

1 1 c,  (?)  (-t)®  + Cp  (?)  t , 11(81) 

where  plus  refers  to  the  liquid  and  minus  to  the  vapor ,5 and-  ::  : t 

C2  (?)  = + ? (c^^'  - 11(82)- 
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Cl  (?)  = 


11(83) 


.(4) 


C (cj^>  - 4'*') 


1 + Cjj  ? (1  - c) 


Q (?.  T = 0)|/P|.  (?) 


where 


Q (?.  ?) 


. <1[RT.(?)]-’ 

^ ^ 

h 


11(84) 


and  is  a parameter  which  characterizes  the  mixing  of  density  and 

41 

composition  difference  across  the  phase  boundary  . The  function 
Q (?,  T)  is  subsequently  transformed  into  a function  of  "c  and  t,  cf.  Eq. 
(12)  of  Ref.  51. 

The  isobutane  parameters  and  are  determined  from  a fit  to 

3 

the  coexisting  density  data  of  Viaxman  and  Gallagher  . The  isopentane 
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parameters  Cj  ^ and  ' are  similarly  obtained  from  the  data  of  Young  . 

Results  are  listed  in  Table  II  4. 

The  liquid  and  vapor  compositions  at  a field  variable  point 

(?.  t)  are  given  by 


X = ? + ? (1  - ?)  H (?,  ?)  . 11(85) 

^ P,  (?) 

where 

H (?.  ?)  = 

Vac 

The  convention  which  defines  x=0  and  x=l  has  been  reversed  from 
Ref.  41;  that  which  is  called  x there  is  called  (1-x)  here.  Equations 
11(84)  - 11(86)  are  analogous  to  Eqs.  11(42),  11(68)  and  11(69)  except 
that  the  term  Q (Tt  ^)/P^  (T)  has  been  added  and  subtracted.  Setting 
p = in  Eq.  11(85)  gives  x^,  and  setting  p = p^  gives  x^. 


)-• 

' T 


d[RT^  (?)]''* 
d? 


fr). 


Q (?■  0) 

P,  (?) 


11(86) 
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In  the  present  model,  H as  a function  of  ^ and  t is  assumed  to  have  ^ 
the  form^^  | 

* I 

_ _ _ d £n  [T  (^] 

H {?.  t = C„  (1  + C,  c)  t . Il  l 

" d? 

I' 

1 

Equations  11(85)  and  11(87)  imply  that  on  the  critical  line,  where  t=0,  ! 

i 

The  critical -line  data  are  fitted  to  polynomials  in  x,  or,  equivalently,' 
C,  according  to  I 


Rf^ 


RT 


1 - X 


RT 


W 


+ x(l-x)  [T^  + (1^2x.)T2  (l"2x-pT-3-}-s  I.;! 


P (x)  xP^^^  (l-x)P^^^  j 

RT^  (x)  ^j(4)  + + ('*“2x)P2  + (1-2x)^P-3]  , 

c c i 

j 

I 

Pc  (x)  = Xp^^^  + (1-x)  + x(l-x)  [p^  + (1-2x)p2  + (1=2x)^  P3],  - - Ilij 

Equation  11(88)  is  analogous  to  Eq.  11(58)  where  T^  ~ ^2~~  "^3  ~ ' • * 

the  present  model,  x is  linear  in  "c  on  the  critical  line  and  the  inverse 
critical  temperature  is  nonlinear  in  "c,  whereas  in  the  model  of  Secs.  II  ' 

3 and  II  4 the^reverse  is  true.  Equations  11(89)  and  11(90)  are  analogous- 
to  Eqs.  11(60)  and  11(59),  respectively.  c ; . : 

Knowledge  of  the  critical  line  and  the  pure  fluid  coexistence-  . - ^ 

properties  determines  all  parameters  of  the  model  except  C^,  and^C^,  -f 
which  are  then  varied  to  optimize  the  fit  of  the  dew-bubble  curves  and -- e:  :: 
mixture  temperature-density  coexistence  curves.  It  turns  outi  however i - 
that  for  some  mixtures  one  or  more  of  these  parameters  may  be  set  equal  to 


zero. 
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_ , 41  51  52 

From  previous  experience  with  this  model  ’ * , Class  I binary 

mixtures  may  be  divided  into  the  following  categories: 

...  34 

(a)  Mixtures  with  narrow  dew-bubble  curves,  e.g.  He  - He,  propane  - 

52  42  51 

butane  , and  most  azeotropic  mixtures  * , for  which  only  a nonzero 

C^,  at  most,  is  necessary; 

41 

(b)  Mixtures  with  medium-width  dew-bubble  curves,  e.g.  nitrogen  - methane 

52 

and  propane  - pentane  , for  which  nonzero  Cm  and  C are  necessary; 

n X 

...  41  . 

(c)  Mixtures  with  wide  dew-bubble  curves,  e.g.  butane  - octane  and 

52 

methane  - ethane  , for  which  nonzero  C^,  C^  and  C^  are  necessary,  and 

(d)  Mixtures  with  extremely  wide  dew-bubble  curves,  e.g.  propane  - 

51  52 

octane  * , for  which  no  set  of  values  of  Cm,  C„,  and  C,  can  be  found 

n x z 

which  provides  an  adequate  fit  of  the  mixture  coexistence  data. 

There  is  substantial  evidence  that  the  mixture  of  isobutane  and 

isopentane  falls  under  category  (a)  above,  i.e.  it  possesses  narrow  dew- 

bubble  curves,  in  contradiction  to  the  results  of  Fig.  II  2.  First, 

dew-bubble  curves  which  are  quite  narrow  have  been  measured  for  the 

54 

mixture  of  normal  butane  and  normal  pentane  by  Kay,  et  al . . Second, 

normal  alkane  mixtures  of  comparable  molecular  mass  ratio,  e.g.  propane  - 
butane,  also  show  narrow  dew-bubble  curves.  Finally,  the  classical 
equation  of  state  for  isobutane  -isopentane,  which  predicts  a critical 
line  approximately  0.1  MPa  too  high  in  pressure  but  which  should  be 
reliable  for  P _<  3.0  MPa,  predicts  narrow  dew-bubble  curves  as  well. 

In  previous  applications  of  this  model,  P-T-p  coexistence  data 
points  were  available,  for  five  or  six  different  compositions,  distributed 
throughout  the  interval  -0.1  < t < 0,  In  the  present  case,  there  are  no 
experimental  data  for  x > 0.51  and  jtj  « 0.1  for  the  data  of  Table  I 7. 
Therefore,  an  alternate  procedure  was  employed.  The  classical  equation 
of  state  was  used  to  generate  pseudo  data  points  at  P = 3 MPa,  and  the 
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1 


critical  line  and  T^(x))  and  the  parameter  Cj^  were  varied  until  a i 

model  was  obtained  consistent  with  both  the  real  and  pseudo  data  points.  I 

I 

Subsequently,  the  classical  equation  was  used  to  generate  further  i 

il 

coexistence  points  for  2.0  MPa  £ P £ 2.8  MPa  in  steps  of  0.2  MPa,  for 

I 

compositions  of  0.2,  0.35,  0.50  and  0.65.  These  results,  and  the  theoretical  " 

j[j 

model  (where  Cj^  = -12),  are  shown  in  Fig.  II  3.  : 

Agreement  is  very  good  for  both  pressure  versus  temperature  and 
temperature  versus  density.  There  is  a small  systematic  trend  toward 
larger  liquid  densities  from  the  model  than  those  predicted  by  the  ; 

classical  equation,  but  the  latter  is  subject  to  refinement  from  additional  i 

high-pressure  thermodynamic  measurements,  and  a discrepancy  of  this 

I 

magnitude  can  easily  be  corrected  by  minor  adjustments  of  the  critical  f 

I: 

line  or  C^,  or  possibly  by  inclusion  of  a nonzero  C^.  The  critical  line  f 

parameters  from  Eqs.  11(88)  - 11(90)  are  listed  in  Table  II  5.  jj 

Also  of  interest  is  the  slope,  dT/dp,  of  the  mixture  temperature  - 
density  coexistence  curves  at  the  critical  line.  In  Table  II  6 we  “ 
compare  dT/dp  as  predicted  by  the  model  with  the  slope  determined  from 
the  data  of  Table  17.  In  the  latter  case  the  slope  is  approximated  by 
the  ratio  AT/Ap  from  the  data  points  farthest  apart  in  density  and 
temperature,  except  that  the  point  x = 0.5073,  T = 437.42  K,  p = 3.629  | 

3 

mol/dm  is  omitted  since  its  temperature  appears  to  be  systematically  high. 

Since  the  experimental  slopes  involve  the  differences  of  very  nearly 
equal  quantities,  agreement  between  theory  and  experiment  is  satisfactory.-- 
We  conclude  that  the  Moldover  - Rainwater  version  of  the  Leung- 
- Griffiths  model,  although  probably  inadequate  for  the  one-phase  region, 
is  capable  of  providing  a satisfactory  fit  in  the  near-critical  region  : 

(-0.1  < t < 0)  for  VLE  of  the  isobutane-isopentane  mixture. 
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Part  III 


Viscosities  of  Hydrocarbons  and  Their  Mixtures 


1.  Experimental  viscosity  data  obtained  under  this  contract 


1.1  Apparatus 

Viscosity  values  are  obtained  from  the  damping  of  a piezoelectric 
quartz  crystal,  5 cm  long,  0.5  cm  in  diameter,  mounted  in  a pressure  cell 
filled  with  the  fluid  of  interest.  The  apparatus  is  described  in  Refs.  12-14. 
The  range  is  100-300  K in  temperature,  at  pressures  up  to  30  MPa.  The 
apparatus  is  most  suited  for  measurements  in  the  saturated  and  compressed 
liquid.  The  viscosity  is  derived  from  the  measured  bandwidth  of  the  resonance 
curve,  Af,  by  means  of  the  relation 


TTf 

M 

Af 

Af 

vac 

P 

S 

f 

f 

vac 

III(l) 


where  p is  the  fluid  density;  M is  the  crystal  mass;  S is  the  crystal  surface 
area  adjusted  for  thermal  expansion  and  hydrostatic  compression;  and  f is  the 
resonant  frequency  of  the  crystal. 

Since  the  damping  is  related  to  the  product  of  viscosity  and  density, 
while  it  is  measured  as  function  of  pressure  and  temperature,  it  is  necessary 
to  have  independent  information  on  the  equations  of  state  of  the  fluids 
studied.  In  the  cases  we  have  studied,  this  information  was  always  available, 
with  an  imprecision  far  smaller  (about  0.2%)  than  that  of  the  measured  bandwidth 
(typically  of  the  order  of  a percent) . The  inaccuracy  of  the  reported 
viscosity  data  is  estimated  to  be  no  larger  than  + 2%,  except  for  the 
butanes,  where  the  imprecision  occasionally  was  as  much  as  + 3%. 
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1.2  Data  for  methane,  ethane  and  propane 


The  viscosity  data  for  methane,  ethane  and  propane  have  all  been 

14-17  * 

published  in  the  archival  literature.  We  display  the  methane  data  in 

Fig.  Ill  1,  those  for  ethane  in  Fig.  Ill  2 and  those  for  propane  in  Fig.  Ill  3. 
These  figures  make  clear  how  the  viscosity  coefficient  of  a typical  non-polar 
fluid  depends  on  density  and  temperature.  It  is  important  to  recall  that 
the  dependence  on  density  is  much  ’simpler’  than  the  dependence  on  pressure, 
especially  at  temperatures  near  the  critical  temperature  of  the  fluid. 

Figure  III  2 shows  that,  at  reduced  densities  smaller  than  about  2.5  (about 
17  mol  L ^ for  ethane) , the  viscosity  typically  increases  weakly  with  both 
temperature  and  density.  However,  at  reduced  densities  larger  than  about 
2.5,  the  viscosity  increases  strongly  with  density  and  decreases  weakly  with 
temperature  at  fixed  density.  We  shall  call  these  density  ranges  the 
’gas-like’  and  ’liquid-like’  density  ranges,  respectively. 

For  the  listings  of  the  original  data  we  refer  to  the  literature  cited. 

1.3  Data  for  nitrogen-methane  and  methane-ethane  mixtures. 

The  viscosities  of  pure  nitrogen,  pure  methane  and  for  approximately 
30-70,  50-50  and  70-30  mol  % mixtures  are  reported  in  Ref.  19.  The  data 
for  the  50-50  mixture  are  displayed  in  Fig.  Ill  4.  The  viscosities  of 
approximately  30-70,  50-50  and  70-30  mol  % mixtures  of  methane  and  ethane 
are  reported  in  Ref.  20.  The  data  for  the  50-50  mixture  are  shown  in  Fig.  Ill  5. 

1.4  Viscosity  data  for  normal  and  isobutane. 

The  viscosity  of  saturated  and  compressed  liquid  normal  butane  was  measured 
at  temperatures  between  140  and  300  K,  and  at  pressures  to  30  MPa.  Densities 

3 

ranged  from  9.8  to  12.8  mol/dm  . The  viscosities  of  the  compressed  liquid 
are  listed  in  Table  III  1,  those  of  the  saturated  liquid  normal  butane  in 
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Table  III  2.  The  data  are  displayed  in  Fig.  Ill  6. 


The  viscosity  of  saturated  and  compressed  liquid  isobutane  was  measured  at 
temperatures  between  120  and  300  K and  at  pressures  to  30  MPa.  Densities 

3 

ranged  from  9.4  to  12.9  mol/dm  . The  viscosities  of  compressed  liquid 
isobutane  are  listed  in  Table  III  3,  those  of  saturated  liquid  isobutane  in 
Table  III  4.  The  data  are  displayed  in  Fig.  Ill  7. 


2.  The  representation  of  the  viscosity  of  hydrocarbon  mixtures. 


2.1  Corresponding-states  model 


A program  for  predicting  the  viscosities  of  hydrocarbon  mixtures  was 

developed  at  NBS  and  is  available  in  the  literature  under  the  acronym  TRAPP. 

This  program  is  based  on  an  extension  of  the  generalized  principle  of 

corresponding  states  discussed  in  Part  I,  Sec.  2.2.  The  reference  fluid  is 

methane.  The  shape  factors  are  not  tailor-made  for  each  substance;  rather, 

they  are  generalized  functions  of  volume  and  temperature  that  depend  on  the 

characteristics  of  the  hydrocarbon  of  interest  through  only  one  parameter, 

28 

the  acentric  factor,  by  a procedure  developed  by  Leach  and  Leland. 

For  the  mixing  rules  used  in  the  program  TRAPP  we  refer  to  the 
literature.  The  adjustable  constants  in  the  combining  rules  are 

set  equal  to  unity. 

The  viscosity  of  the  mixture  at  specified  pressure,  temperature  and 
composition  is  obtained  in  the  following  steps.  Starting  points  are  the 
equation  of  state  and  viscosity  formulations,  in  terms  of  density  and 
temperature,  for  the  reference  fluid.  The  program  begins  by  finding  the 
density  of  the  mixture  and  the  equivalent  density  and  temperature  of  the 
reference  fluid,  using  the  generalized  principle  of  corresponding  states 
for  the  mixture.  Once  the  densities  are  known,  a simple  expression  connects 
the  viscosity  of  the  mixture  to  that  of  the  reference  fluid  at  the  corresponding 
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density  and  temperature.  The  program  TRAPP  claims  an  accuracy  of  about 
seven  percent  in  the  prediction  of  transport  properties  of  hydrocarbon 
mixtures . 

2.2  Comparison  of  the  corresponding- states  model  with  experimental  data 
for  methane,  ethane  and  propane 

Figure  III  8 compares  our  data  for  methane  at  100  K with  a correlating 
equation  for  methane  only.^^  Figures  III  9 and  10  compare  our  data  for 
saturated  and  compresed  liquid  ethane  and  propane,  respectively,  with 
the  corresponding-states  model.  The  densities  shorn  correspond  to 
reduced  densities  between  1.6  and  3*2.  The  comparisons,  along  with 
others  in  the  publications  cited,  support  the  observation  that,  in  the  . _ 
’gas-like’  density  range,  most  of  the  differences  between  our  measurements 
and  the  model  are  smaller  than  7 percent,  and  that  the  model  is  likely 
to  be  satisfactory  for  many  engineering  calculations.  However,  at 
reduced  densities  larger  than  about  2.5,  we  find  that  in  most  cases  the 
differences  are  substantially  larger  than  7 percent,  and  that  they 
typically  increase  strongly  with  density. 

2.3  Comparison  of  the  corresponding- states  model  with  data  for  the  mixtures 
of  nitrogen  and  methane,  and  those  of  methane  and  ethane 

The  viscosity  measurements  on  pure  methane,  ethane  and  propane  have 

revealed  some  of  the  strengths  and  weaknesses  of  the  corresponding- 

states  model  applied  to  pure  fluids.  The  viscosity  measurements  that  we 

19  20 

have  performed  on  nitrogen-methane  and  methane-ethane  mixtures  ' ' 

serve  to  answer  the  question  of  whether  the  model  degrades  if  applied  to 
mixtures . 

The  case  of  the  nitrogen-methane  mixtures  is  not  clear-cut.  As  is 
shown  in  Fig.  Ill'  11,  the  measured  viscosities  of  the  30-70,  50-50  and  70-30 

Sfi  


mol  % mixtures  depart  systematically  from  the  model,  but  the  departures  are 
generally  within  5%.  The  concentration-dependence  at  fixed  molar  density 
and  temperature  is  not  particularly  well-predicted  (Fig.  Ill  12)  but  the 
problem  is  caused  by  poor  prediction  of  the  pure-nitrogen  viscosity,  a 
result  perhaps  not  surprising  given  that  the  model  was  developed  for  use 
in  hydrocarbon  mixtures. 

The  methane-ethane  mixtures  form  a more  interesting  test  case.  In 
Fig.  Ill  13,  we  show  a comparison  of  pure-methane  and  pure-ethane  viscosities 
with  the  corresponding-states  model.  Up  to  15%  departures,  of  opposite 
signs  for  the  two  fluids,  are  found  at  the  highest  densities.  Note  that  the 
reference  viscosity  values  of  the  model  have  been  shifted  with  respect  to  the 
correlation  for  pure  methane  (Fig.  Ill  8)  in  order  to  represent  viscosities 
of  the  heavier  hydrocarbons  more  closely.  In  Fig.  Ill  14,  we  shown  viscosities 
of  the  approximately  35-65  mol  % methane-ethane  mixture  over  the  entire  range  of 
densities.  Since  at  the  lowest  densities  the  agreement  is  within  5%,  and  the 
largest  departures  are  within  15%,  there  seems  in  this  case  to  be  no  further 
degradation  of  the  model  if  extended  to  mixtures. 

2.4  Comparison  of  the  corresponding-states  model  with  data  for  normal  and 

isobutane 

Of  most  interest  to  the  objectives  of  this  project  are,  of  course,  our 
data  on  the  viscosities  of  normal  and  isobutane. 

Measured  normal  butane  viscosities  were  compared  with  the  generalized 

corresponding  states  model.  Measured  and  calculated  saturated  liquid 

viscosities  are  shoTcn  in  Fig.  Ill  15.  The  differences  between  the  measured 

and  calculated  values,  for  both  compressed  and  saturated  liquid,  are  shown  in 

Fig.  Ill  16.  At  lower  densities  (<  2.5  p ) the  differences  are  smaller  than 

— c 

5 percent,  which  is  typical  of  other  fluids  we  have  examined  in  this  reduced 
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I 

I 

! 

density  range.  At  densities  larger  than  2.3  P the  differences  increase  with  : 

i 

density  to  a maximum  difference  of  about  25  percent,  which  is  also  typical  of  j 

'i 

our  comparisons  for  other  fluids  in  this  reduced  density  range.  j' 

Comparisons  of  measured  saturated  liquid  isobutane  viscosities  with 

I ; 

the  same  extended  corresponding  states  model  are  shown  in  Figs.  Ill  15,  l6,  .j 

Differences  between  measured  and  calculated  viscosities,  for  both  saturated  !' 

i) 

and  compressed  liquid  data,  are  shown.  The  differences  increase  with 
density  from  about  15  percent  at  the  lowest  densities  to  over  100  percent 
at  the  highest  densities.  At  the  highest  densities  the  differences  are  ; 

much  larger  for  isobutane  than  for  normal  butane,  which  is  somewhat  1 

surprising  considering  that  they  are  isomers.  I' 

Our  measurements  under  this  contract  have  therefore  revealed  an  isomer  li. 

I 

effect  on  viscosity  that  is  much  larger  than  the  generalised  corresponding-states  jj 
model  is  capable  of  giving.  This  suggests  that  the  program  TRAPP  should  be  j 

used  with  caution  at  densities  higher  than  2p  for  the  isobutane-isopentane  j 

mixtures  of  interest  in  geothermal  power  cycles.  j 

I 

2.5  An  improved  formulation  of  hydrocarbon  viscosities 

Hildebrand^^  has  shown  that  the  equation, 

(J)  = B(V  - , IlKl) 

where  (|)  = n ^ is  called  the  fluidity,  B is  an  empirical  parameter,  and  is 

the  estimated  molar  volume  at  (J)  = 0,  can  correlate  a large  amount  of  viscosity 

data  for  liquids  at  low  reduced  temperatures.  It  is  known  that  equation  III(l) 

is  a satisfactory  approximation  only  within  limited  ranges  of  ten^Deratures 
57 

and  pressures.  Nevertheless,  it  provides  a relatively  simple  way  to 
characterize  the  viscosity  behavior  of  liquids  within  these  ranges. 

Eq.  III(l)  suggests  that  a plot  of  the  inverse  viscosity  versus  the  volume 
should  be  linear.  The  slope  of  such  a plot  should  give  the  parameter  B,  the 
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extrapolated  intercept  with  the  r|  =0  axis  the  parameter  V^.  We  have  made 
such  fluidity  plots  for  the  liquids  ethane  (Fig.  Ill  17),  propane  (Fig.  Ill  18), 
and  normal  and  isobutane  (Fig.  Ill  19).  It  is  obvious  that  these  representations 
are  simpler  and  more  revealing  than  the  viscosity  vs.  density  plots.  In  the 
cases  of  ethane  and  propane,  the  constant  seems  to  be  virtually  independent 
of  temperature,  while  the  saturation  boundary  is  closely  approximated  by  a 
straight  line. 

The  dependence  of  the  fluidities  of  the  butanes  on  molar  volume  is 

shown  in  Fig.  Ill  19.  The  fluidities  of  the  compressed  liquid  increase 

linearly  with  volume  at  fixed  temperature;  the  fluidities  of  the  saturated 

liquid  also  increase  linearly  with  volume  except  at  the  lowest  temperatures. 

The  parameter  B is  clearly  dependent  on  temperature  at  low  reduced 

temperatures.  Estimates  of  V^,  obtained  by  linear  extrapolation  of  the 

isotherms  to  (})  = 0,  are  also  slightly  dependent  on  temperature.  The 
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literature  indicates  that  the  compressed  liquid  isotherms  would  probably 
be  non-linear  at  volumes  smaller  than  we  have  measured.  The  behavior  shown 
in  Fig.  Ill  19  suggests  that  equation  III(l),  if  modified  to  include 
dependences  on  temperature,  can  be  used  for  correlating  our  data.  We  have 
therefore  examined  our  data  making  use  of  equation  III(l),  modified  to 
include  dependences  of  B and  on  temperature.  For  normal  butane,  the 
temperature  dependences  are  described  by 

B(T)  = 0.287  exp  (-0.198  x 10^/T^)  III(2) 


and 


V (T)  = 0.074  + 3.33  x 10  ^ T , 
o 

where  T is  in  kelvins,  in  L mol  ^ and  B in  jiPa  ^‘s  (L*mol  ^)  , 

shows  a comparison  of  our  data  with  equations  III(l),  (2)  and  (3). 


III(3) 
Figure  III  20 
Most  of  the 


differences  between  measured  and  calculated  viscosities  are  smaller  than  3 percent. 


For  isobutane  the  temperature  dependences  are  described  by 

B(T)  = 0.420  exp  (-0.354  x 10^/T^) 

and 

V (T)  - 0.0698  X 10"^T  . III(5) 

o 

Figure  III  21  shows  a comparison  of  our  data  for  isobutane  with  equations 
III(l),  (4)  and  (5).  As  for  normal  butane,  most  of  the  differences  between 
measured  and  calculated  viscosities  are  smaller  than  3 percent. 

We  conclude  that  a simple  Hildebrand  equation  with  few  adjustable  parameters 
seems  to  be  effective  in  describing  hydrocarbon  viscosities  in  the  range  from 
2 to  3 p^.  However,  unlike  TRAPP,  the  Hildebrand  equation  is  not  global  and 
not  predictive. 


This  identification  does  not  imply  recommendation  or  endorsement  by  the 
National  Bureau  of  Standards,  nor  does  it  imply  that  the  materials  or 
equipment  identified  are  necessarily  the  best  available  for  the  purpose. 
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Scaled  equations  for  one-component  fluids 
B.l  Reduced  thermodynamic  quantities 
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B.3  Fundamental  equations 
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B.4  Derived  thennodynamic  quantities 
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B,5  Parametric  equations 
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B.6  Thermodynamic  properties  in  one-phase  region 
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B.7  Thermodynamic  properties  in  two-phase  region 
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Scaled  equations  for  fluid  mixtures  (Sec.  3 of  Part  II) 
C,1  Reduced  thermodynamic  quantities 
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C.3  Parametric  representation 


T = Jl(^)T  = r(l-b^e^)  (C.17)  I 

i 


h = r^‘^ae(l-e^) 


(C.18) 


AP  = r^^^'''^^ak(pQ  + + p^9^)  (C.19) 


C4 


■i  + 

‘i  C.4  Derived  chermodynamic  quantities 

i 

(C.20) 

(C.21) 

(C.22) 

Mole  fraction  of  isopentane: 

Pc  Q 

X = — = ^ + 5(1-5)  — (C.23) 

P P 


Densities 


“ P4  ' '’h 


. (1-Op  - ?(l-OQ 


= 5P  + ?(1-4)Q 


Notation 


C5 


Specific  energy: 


(Co24) 


Specific  enthalpy: 


H = U 


(C.25) 


Isochoric  thermal  pressure  coefficient: 


T 

c 


(P 


+ 


(-Co  26) 


Compressibility : 


1 

P 


(Co27) 


C6 


Response  function 


( 


3x 

3A 


3x 


1 T 

7(4)  r 

P 


(C.28) 


Specific  heat  at  constant  volume: 


C 

vx 


P 


(4)^  (4) 
c c 


E 

P F 


(C.29) 


Specific  heat  at  constant  pressure: 


C 

px 


C 

vx 


(C.30) 


C7 


Auxiliary  functions: 


Q=P-PT  -p(H-HT) 
? T cc  ^ T c? 


(C.31) 


D = X (P  p -pp  )-x  (P,  P , - P P,  . ) + X (P,  P , - P P,  , . (C.32) 

T 5 Th^  T h 5 hh  C h xh  x hh)  ^ ^ 


E = X (P  P - P P . ) - X (P  , P , - P P,  , ) + X (P\  - P P,,  ) 

h XX  5x  xh  X xh  ^x  hh  i,  xh  xx  hh 


^ ^«5T\h  - ^«tT"hh] 


(C.33) 


F = [1  + (1-20  ^ + 


hh 

P P 


I pK^  f-H  + 2H  T + H T - H 1 

^ hh  '■  CT  cc  T c??  XX  c^-' 


+ 


P.  . [P  - 2P  T - P T + P 1 
hh  TC  T XX  cc-* 


cc 


[P  , - P 1^1 

(;h  xh  c?-*  ^ 


(C.34) 


g = p + {l~20  Q + 


?(1-?)  {-p(H  - 2H  f + H - H T ) 

CT  CC  XX  ct  X cCC^ 


44  4T  c?  XX  cc  X cC^ 


+ (P  - 2P  T + P T - P T ) 
CT  c?  XX  C4  T c^C 


(P  - P T )(P  ^ - P ) 

Z,  X cc  4h  xh 


1 ^ -V  ^ ,v  2 

P^^,(P  - P T ) 

-2  hh  i;  X c4 


(C.35) 
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Appendix  D 


Parameters  for  isobutane  - isopentane  mixtures  near  the  critical  line 
for  the  scaled  model  of  Part  II,  Sections  3 and  4 


"cO  • ^cl  ' 


cO 


p , » -0.15842 

cl 


p^2  = 0.06 


^00  " ^ 


Pqj  » -0.17735 


Pq2  = 0.16 


03 


P^Q  « 5.8840  Pj^  » -0.15089 


Pl2  = 7 


P20  = -15.789  P2J  = -3.1202 


^30  “ ° 


P^^  = 23.089 


^00  “ ^ 


= -0.022491 


- 5.8840  = 0.92834 


H2Q  » -28.659  H2j  = -23.22 


H30  * -25.543  = 3.7432 


* 1 q,  ■ -0.15401 


JIq  = 1 ^1  “ 0.11264 


D1 


-0.03 


I ' 


k}  t.  ft'rroi  ' drf:*:.  xolv '..: 


: ri'  ;,'■  I ^ 

■'r 


■>.,  ^>  - , . . -■£ 
i >’w 

' ' .. . • f -'^^1 

Bi?  l .0*1*'  ' '*•  "‘ij  \ ' I * .,  Q "':  j 

-■  ■ '5  ■ rX'2  -'-  •\' '■ -ft  0.a^--  ■'  ■'- 
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Appendix  E 


Tables  of  thermodynamic  properties  of  near-critical  isobutane-isopentane  at  x=0.1 

In  this  Appendix  we  present  tables  of  thermodynamic  properties  for 

isobutane  - isopentane  mixtures  calculated  from  the  scaled  fundamental 

equation.  Specifically  we  list  the  pressure  in  MPa,  the  thermal  pressure 

coefficient  (3P/3T)  in  Mpa/K,  the  isothermal  compressibility 

P 9^ 

= p ^(3p/3P)^  ^ in  1/MPa,  the  specific  heat  in  kJ/kg  K,  the  specific 
heat  Cp  in  kJ/kg  K and  the  specific  enthalpy  H in  kJ/kg  as  a function  of  the 
density  p in  kg/m  and  temperature  in  K for  two  compositions,  namely  x =•  0.1 
and  X = 0.5,  where  x represents  the  mole  fraction  of  isobutane. 
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Table 


Table 


-Table 
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Table 


Pressure  in  MPa  as  a function  of  p in  kg/m  and  T in  K 
for  X = 0.1  3 


E2.  Thermal  pressure  coefficient  (3p/3T)  in  MPa/K  as  a function 

« P 

of  p in  kg/m"^  and  T in  K for  x = 0.1 

E3,  Isothermal  compressibility  = p ^(8p/3P)^  ^ in  1/Mpa  as  a 

3 * 

function  of  p in  kg/m  and  T in  K for  x * 0«1 

3 

E4,  Specific  heat  C-,  „ in  kJ/kg  K as  a function  of-p  in  kg/m  and 

V 

T in  K for  x = 0.1 

3 

_ _ in  kJ/kg  K as  a function  of  p in  kg/m  and 

E5.  Specific  heat 

T in  K for  x - 0. 1 

3 

E6.  Specific  enthalpy  H in  kJ/kg  as  a function  of  p in  kg/m  and 
T in  K for  x - 0.1 
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E2.  Thermodynamic  properties  of  isobutane  - isopentane  at  x = 0,5 


Table  E7.  Pressure  in  MPa  as  a function  of  p in  kg/m^  and  T in  K 
for  X = 0.5 


Table  E8,  Thermal  pressure  coefficient  (3p/3T)  in  Mpa/K  as  a function 

_ P 

of  p in  kg/m"^  and  T in  K for  x = 0.5 


Table  E9. 


-1 


Isothermal  compressibility  Kj.  = p (3p/3P)^  ^ in  1/MPa  as  a 

■3  * 

function  of  p in  kg/m  and  T in  K for  x = 0,5 


Table  ElO,  Specific  heat  C in  kJ/kg  K as  a function  of  p in  kg/m^  and 

V 

T in  K for  x = 0,5 


Table  Ell.  Specific  heat  C in  kJ/kg  K as  a function  of  p in  kg/m"^  and 

p ,x 

T in  K for  x = 0.5 


Table  E12. 


Specific  enthalpy  H in  kJ/kg  as  a function  of  p in  kg/m  and 
T in  K for  x = 0.5 
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Table  I 1 


Supercritical  PVT  data  for  isobutane 


T = 

423.15  K 

T = 

448. 15K 

P,  MPa 

p,  kg/m^ 

P,  MPa 

P,  kg/ 

20.750 

457.9 

10.024 

349. 

4.754 

257.2 

5.770 

196. 

10.840 

406.9 

4.377 

110. 

4.600 

228.5 

3.018 

61. 

3.977 

128.3 

10.024 

349. 

6.426 

349.3 

5.770 

196. 

4.452 

196.2 

4.377 

110. 

3.747 

110.2 

6.757 

261. 

5.403 

314.1 

5.060 

146. 

4.352 

176.4 

3.675 

82. 

3.571 

99.91 

7.480 

292. 

15.982 

438.4 

5.328 

164. 

4.690 

246.3 

3.946 

92. 

4.078 

138.3 

uncertainty  0.01%  0. 


6 

3 

3 

97 

6 

4 

3 

2 

7 

39 

5 

2 

24 

1% 


T1 


Table  I 2 


Virial  coefficients  of  isobutane 


lO^B. 

lO^C, 

lO^D. 

lO^^E. 

P 

max 

T,K 

m^/kg 

(m^/kg)^ 

(ra^/kg)^ 

(m^/kg)"^ 

kg/m 

377.594 

-647.4  + 0.7 

9,86  + 0.16 

41 

394.261 

-592.2  + 0.2 

10.66  + 0.08 

81 

423.15 

-501,5  + 0.3 

8.91  + 0.07 

7.9  +0.8 

-3.0  + 0.3 

150 

448.15 

-440,3  + 0.1 

8.62  + 0.01 

100 

B,  second  virial;  C,  third  virial;  D,  fourth  virial,  E,  fifth  virial» 
Virial  coefficients  correspond  to  the  coefficients  of  a finite  density 
series  representation  of  the  compressibility  factor. 


T2 


Table  I 3 


Vapor  pressure  of  isobutane 


T,K 

P ,MPa 

298.15 

0.3500 

303.15 

0.4043 

308.15 

0.4641 

313.15 

0.5304 

318.15 

0.6032 

323.15 

0.6836 

328.15 

0.7725 

333.15 

0.8683 

Uncertainty 


T,K 

P,MPa 

333.15 

0.8676 

343.15 

1.0867 

353.15 

1.3432 

363.15 

1.6408 

373.15 

1.9855 

383.15 

2.3819 

393.15 

2.8361 

398.15 

3.0879 

0.002 

0.01% 

T3 


Table  I A 


Isochoric  PpT  data  near  the  critical  line 
65.07/34.93  mol  % isobutane/isopentane  mixture 


T,K 

P,MPa 

T,K 

P,MPa 

p = 251. 

9 kg/m^ 

p - 241 

3 

.2  kg/m 

157.535 

3.801 

156.582 

3.729 

157.424 

3.795 

156.520 

3.726 

157.362 

3.791 

156.444 

3.721 

156.558 

3.804 

156.468 

3.722 

157.415 

3.796 

156.328 

3.714 

157.918 

3.766 

156.113 

3.  701 

156.635 

3.747 

155.708 

3.676 

156.563 

3.742 

155.536 

3.666 

156.474 

3.735 

p = 230. 

8 kg/m^ 

156.389 

3.730 

156.333 

3.  726 

156.503 

3.703 

156.330 

3.726 

156.369 

3.696 

156.286 

3.691 

155.811 

3.693 

155.974 

3.675 

155.761 

3.693 

155.865 

3.668 

155.697 

3.690 

155.719 

3.692 

p = 216. 

7 kg/m 

155.628 

3.686 

156.606 

3.700 

155.497 

3.678 

156.533 

3.696 

155.425 

3.673 

155.972 

3.665 

155.245 

3.661 

156.233 

3.681 

uncertainty 

156.364 

3.689 

0.02 

0.004 

156.413 

3.691 

156.458 

3.693 

156.529 

3.698 

T4 


Table  I 5 


Critical-line  data  in  isobutane/isopentane  mixtures 


X,  mol  % iCf, 

T 

p 

. P 

5 

c 

c 

c 

K 

3 

mol/ dm 

MPa 

0 

407.84 

3.880 

3.629 

0.2001 

419.92 

3.781 

3.677 

0.3493 

428.71 

3.662 

3.650 

0.5073 

437.35 

3.649 

3.615 

1 

460.51 

3.265 

3.371 

uncertainty 

0.02 

1% 

0.004 

T5 


Table  i 6 


Vapor  pressure  of  isopentane 


T,K 

P ,MPa 

exp 

400.529 

1.2477 

411.407 

1.5208 

411.390 

1.5222 

425.513 

1.9387 

437.593 

2.3593 

447.995 

2.7761 

454.576 

3.0728 

458.582 

3.2698 

459.562 

3.3192 

460.504 

3.3698 

0.02 

0.004 

Fitting  equation  11(25). 


P , , MPa  P -P 
calc  exp 


1.2467 

.0010 

1.5230 

-.0022 

1.5226 

-.0004 

1.9373 

.0014 

2.3570 

.0023 

2.7770 

-.0009 

3.0755 

-.0027 

3.2718 

-.0020 

3.3218 

-.0026 

3.3707 

-.0009 

uncertainty 


(P  -P  , )/P  , % 

exp  calc  exp 

.08 

-.15 

-.03 

.07 

.10 

-.03 

-.09 

-.06 

-.08 

-.03 


T6 


Table  I 7 


Dew-bubble  curve  data 


X 

T,K 

3 

p,  mol/dm 

P , MPa 

mol  % iC^ 

0.2001 

419.75 

3.993 

- 

0.2001 

419.88 

3.832 

- 

0.2001 

420.01 

3.673 

3.676 

0.3493 

427.99 

4 

4.336 

- 

0.3493 

428.14 

4.161^ 

- 

0.3493 

428.41 

3.996^ 

3.661 

0.3493 

428.56 

8 

3.826 

3.658 

0.3493 

428.71 

9 

3.661 

3.650 

0.3493 

428.89 

3.438^ 

3.650 

0.5073 

436.88 

4.005^ 

- 

0.5073 

437.11 

3.825^ 

3.62^ 

0.5073 

437.28 

3.662^ 

3.637 

0.5073 

437.42 

3.629^ 

3.62 

uncertainty 

0.001 

0.02 

1% 

0.004 

(x  = 0.20  and  0.35) 
0.01  (x  = 0.50) 


T7 


Table  I 8 


VLE  data 


T 

X,  mol  % iC 

288.890 

0,4968 

0.5010 

0.5097 

298.337 

0,4968 

0.5010 

308.327 

0.4968 

0.5010 

318,496 

0.4968 

0.5010 

0.5097 

328.480 

0,4968 

0.5010 

P,MPa 

p,  mol /dm' 

0.1685 

9.3110 

0.1653 

9.3041 

0.1632 

9.2790 

0.2270 

9.1483 

0.2238 

9.1449 

0.3045 

8.9542 

0.2994 

8.9469 

0.4029 

8.8137 

0.3965 

8.7712 

0.3894 

8.8082 

0.5145 

8.6170 

0.5139 

8.5955 

T8 


Table  I 9 


Parameters  of  the  dimensionless 
Helmholtz  function  of  isobutane 

Primary  reference  constants 
T*  = 407. 84K 
p*  = 3.8796  mol/dm^ 

P*  = 3.6290  MPa 

Secondary  reference  constants 
A**  - = 0.93541  kJ/mol 

p* 

S**  - = 2.29356  J/(mol.K) 

P* 

a**  - fe)  = 126.86  m/s 

p« 

Physical  constants 
M = 0.0581243  kg/mol 
R = 8.31441  J/(mol.K) 


T9 


Table  I 9 (Continued) 


Coefficients  of  dimensionless  Helmholtz  function 


XO  = 7.9688113 


AGO  = 27.849029 
AOl  = 58.138236 
A02  = -.050937 
A03  = 2.525250 
A04  = 346.65236 


A05  = 2.5515510 

A06  = -53.012165 
A07  = 3.343460 

A08  = -.119451 

A09  = -38,124442 


YO  = 

.15388314 

A20 

= 

.30020353+001 

Y1  = 

■.39169870-001 

A21 

-.61529971+001 

Y2  = 

■.25198404-003 

A22 

-.14570002-001 

Y3  - 

.98801205-006 

A23 

= 

.13342155 

A24 

-.90043710-004 

ZO  = 

' • 

38796166 

AlO 

= 

3.6250 

Bll 

= 

-.96153074+001 

B23 

-.50009149+002 

B21 

= 

.27935713+002 

B43 

.23153999+003 

B41 

= 

-.12569635+003 

B63 

= 

-.38080769+003 

B51 

= 

.54406550+003 

B83 

.26120687+003 

B61 

= 

-.47948565+003 

B14 

= 

-.22934154+002 

B81 

= 

.14134133+003 

B64 

= 

-.14503027+002 

B12 

= 

-.12372626+002 

B15 

= 

-.10167777+002 

B32 

= 

-.34731447+002 

B25 

= 

.30142576+002 

B52 

= 

-.57569010+003 

B55 

= 

-.33549797+002 

B62 

= 

.53210066+003 

B85 

= 

.25502886+002 

B72 

= 

.41502454+003 

B26 

= 

-.53441617 

B82 

= 

-.42359614+003 

B86 

= 

.37213690-001 

B13 

= 

.58118955+002 

TIO 


Table  I 10 


Coefficients  determining  shape  factors  9 and  (p 


9 

<P 

^1 

1 

1 

^2 

+0.0065 

*>2 

-0.021 

^3 

-0.015 

^3 

-0.065 

Til 


Table  I 11 

Critical  parameters  of  isobutane  and  isopentane 


T,^  = 407.84  K 
4 

= 0.2578  dm^/mol 
4 

= 3.629  MPa 
4 

Z,^  = 0.2759 
4 


= 460.51  K 

= 0.3079  dm  /mol 

P.^  = 3.3707  MPa 
5 

= 0.2711 


T12 


Table  I 12 

Combining-rule  parameters 


k - 0.995 

I = 1.002 


and  % 


T13 


Table  I 13 


Parameters  of  the  revised  and  extended 
scaled  thermodynamic  potential  of  near-critical  isobutane 

Critical  Parameters 

T = 407.84  K P = 3.6290  MPa 

c c 

= 225.5  kg/m^ 

Critical  Exponents 
8 = 0.325  = 0.50 

6 = 4.82 

Parameters  in  Scaling  Functions 

a = 22.0163  c - -0.0096833 

k = 1.19385  b^  = 1.3757 

o 

k^  = 0.50552 

Pressure  Background  Parameters 

P^  = 5.8858  P^^  = -0.068209^ 

P^  = -22.0805 

Thermal  Background  Parameters 

y = -4.9535  y„  = -32.2295 

c 2 

y^  = -21.6912  y^  = -33.  5271 


T14 


Table  II  1 


Par^lmeters  for  the  simple  scaled  equtions  of  isobutane  and 

(Sec.  2 of  Part  II) 

isobuCane 


Molar  mass 

M 

sv 

0.0581243  kg/mol 

M 

Critical  parameters 

^.(4) 

c 

3 

407.84  K 

^(5) 

c 

s 

3.8796  mol/dm^ 

p(4) 

c 

s 

3.62906  X 10^  Pa 

p(5) 

c 

Parameters  in  scaling 

6 

= 

0.355 

3 

functions 

5 

4.352 

6 

a 

3 

20.7234 

a 

k 

= 

1.44675 

k 

b2 

3 

1.3909 

b2 

Pressure  background 

p(4) 

1 

3 

5.8840 

-^(5) 

1 

parameters 

2 

3 

-15.7886 

-p(5) 

-(4) 

3 

0 

^p(5) 

3 

3 

Caloric  background 

“c 

3 

1 

“c 

parameters 

-(4) 

3 

•5(4) 

1 

-(5) 

-(4) 

X 

-28.6595 

;(4) 

3 

-25.5431 

Range  of  validity 

407. 

6 

K < T < 420.1  K 

460 

isopentane 

isopentane 

0.0721512  kg/mol 

460.51  K 
3.265  mol/dm^ 
3.371  X 10^  Pa 

0.355 

4.352 

20.7234 

1.44675 

1.3909 

6.1720 

-18.0281 


19.4960 

1 


-41.6268 

-15.4912 

K < T < 474.3  K 


Table  II  2 


Calculated 

critical-point 

parameters  of  isobutane 

-isopentane 

mixtures  as 

; a function  of 

the  mole  fraction  x of 

isopentane 

(Sec.  3 of 

Part  II). 

X 

T 

Pc 

Pc 

K 

mol/dm^ 

MPa 

0 

407.84 

3.880 

3.629 

0o2001 

419.40 

3.774 

3,677 

0e3493 

429.14 

3,670 

3.662 

0e5 

437,71 

3.569 

3,616 

Qe5073 

438.09 

3,564 

3,613 

1 

460.51 

3.265 

3.371 

Table  II  3 


Note : 


Differences  between  thermodynamic  property  values  Pg^obal 
calculated  from  the  global  equation  (Sec.  2.2  of  Part  I) 
and  thermodynamic  property  values  calculated  from 

the  scaled  fundamental  equation  (See.  3 of  Part  II)  for  an 
equimolar  mixture  of  isobutane  and  isopentane  (x  = 0.5)  at 
T = 452  K. 


Quantity  listed: 


^^global  ^scaled^^^scaled 


P 


F » P 


3 

mol/ dm 

% 

% 

% 

2.702 

6.1 

0.8 

- 5.6 

2.886 

5.8 

3.4 

-10.9 

3.070 

5.5 

5.8 

-15.1 

3.255 

5.1 

8.0 

-17.9 

3.439 

4.6 

10.0 

-19.5 

3.623 

4.2 

11.8 

oo 

. 

1 

3.807 

3.7 

13.2 

-19.3 

3.992 

3.2 

14.3 

-18.1 

4.176 

2.7 

15.1 

-16.5 

4.360 

2.2 

15.5 

-14.8 

= 437.71  K, 

3 

p ® 3.569  mol/dm  , M = 

0.065137  kg/raol. 

Table  II  4 


Coefficients 


and 


(Sec.  5 of  Part  II) 


lo984 

-0.874 

30.22 

5.874 

-25.45 

-3.287 


2.001 

-0.800 

30.22 

6.249 

-20.94 

16.69 


T18 


Table  II  5 


Critical  line  parameters  (Sec.  5 of  Part  II) 


P. 

1 


Pi 


-1.120  X lO"^  1.023  X lO"^ 


2.875  X 10 


-2 


■1.624  X 10“^  1.196  X lO"^ 


3.490  X 10 


-2 


3.808  X 10"^  5.202  x lo"^ 


9.470  X 10 


-2 


Table  II  6 


Critical 


X 

0.2 


0.35 


0.5073 


slope  of  temperature-density  coexistence  curves  (Sec.  5 of  Part  II) 


, K*dm^/mol  (expt.) 


■0.81 


•1.00 


■1.17 


, K*dni^/mol  (model) 


-0.70 


-0.96 


-0.93 


T20 


Table  III  1 


Viscosities  of 

compressed  liquid  normal 

butane 

P,  hPa 

p,  mol#  L“^ 

T = 300.00  K 

33.9166 

10.548 

221.4 

30.6994 

10.495 

217.9 

28.1542 

10.450 

212.9 

10.389 

205.9 

21.5095 

10.327 

201.3 

17.7614 

10.252 

192.3 

14.4315 

10.181 

188.2 

10.2694 

10.086 

175.8 

7.6515 

10.021 

170.8 

4.8175 

9.947 

164.3 

1.7289 

9.859 

157.8 

T = 250.00  K 


30.0269 

11.201 

337.3 

27.0472 

11.162 

327.3 

24.6018 

11.130 

322.7 

21.4561 

11.087 

316.2 

17.9370 

11.037 

305.9 

14.6209 

10.988 

299.3 

11.7832 

10.945 

292.5 

7.9407 

10.884 

283.6 

5.1172 

10.837 

275.9 

1.8463 

10.781 

267.1 

T21 


Table  III  1 (Continued). 


P,  MPa 

p, 

T 

= 200.00 

K 

31.0293 

11.928 

646.1 

30.2055 

11.920 

632.7 

27.9013 

11.898 

618.5 

24.6365 

11.865 

606.0 

21.3166 

11.831 

587.0 

17.9753 

11.797 

573.8 

14.4286 

11.759 

560.9 

11.1406 

11.723 

550.1 

7.9864 

11.687 

533.8 

5.1249 

11.654 

519.1 

2.0330 

11.618 

507.6 

T 

II 

»— * 

00 

o 

e 

o 

o 

K 

31.2351 

12.219 

880.7 

28.0523 

12.191 

851 .9 

24.7734 

12.162 

835.6 

21.1814 

12.129 

813.8 

17.7851 

12.098 

797.9 

14.3967 

12.065 

782.3 

11.3084 

12.035 

760.8 

8,1990 

12.004 

745.4 

5.4201 

11.976 

729.9 

2.1837 

11.942 

703.7 

T22 


Table  III  1 (Continued) 


P,  I'lPd  p,  inol»L“^  jUPd^S 

T = 160.00  K 


32.2269 

12.519 

1381.8 

30.0563 

12.503 

1368.3 

27.0460 

12.479 

1317.8 

24.4627 

12.458 

1282.6 

21.1290 

12.430 

1248.5 

16.9623 

12.395 

1230.0 

14.4448 

12.374 

1196.0 

11.4659 

12.348 

1169.8 

8.0570 

12.317 

1130.3 

4.6655 

12.286 

1098.3 

1 .4602 

12.257 

1073.4 

T = 150.00  K 


29.1128 

12.644 

1809.5 

26.6361 

12.625 

1785.9 

24.0565 

12.605 

1736.8 

21.4113 

12.585 

1663.2 

17.8222 

12.556 

1641.7 

14.6675 

12.530 

1562.4 

7.9562 

12.474 

1512.4 

5.0882 

12.450 

1443.8 

1.7692 

12.420 

1400.2 

Table  III  1 (Continued) 


P,  MPa  p,  mol'L"^  » uPa.s 


T = 140.00  K 


30.2935 

12.802 

2565.4 

27.7097 

12.784 

2486.9 

24.6230 

12.761 

2440,3 

21.2652 

12.736 

2384.6 

17.4590 

12.707 

2320.1 

15.4894 

12.692 

2264.0 

12.8373 

12.672 

2199.7 

10.2410 

12.651 

2117.4 

7.4174 

12.628 

2055.3 

4.6766 

12.606 

2012.5 

1.5161 

12.579 

1952.4 

T24 


Table  III  2 


Viscosities 

of  saturated 

liquid  normal 

butane 

T,  K 

p,  inol»  L' 

•1 

/iPa.s 

136.00 

12.634 

2234.9 

138.00 

12.601 

2013.1 

140.00 

12.566 

1899.7 

145.00 

12.486 

1592.1 

150.00 

12.405 

1332.7 

155.00 

12.324 

1200.4 

160.00 

12.243 

1047.2 

170.00 

12.081 

816.9 

180.00 

11.920 

690.3 

190.00 

11.757 

571.5 

200.00 

11.593 

485.2 

220.00 

11.262 

368.2 

240.00 

10.922 

288.7 

250.00 

10.749 

259.8 

260.00 

10.572 

233.4 

280.00 

10.204 

189.0 

t 

l| 


T25 


Table  III  3 


Viscosities  of  compressed  liquid  isobutane 


P,  MPa  p, 


T = 300.00  K 


32.0838 

10.250 

223.8 

28.5975 

10.182 

215.6 

25.1765 

10.112 

206.3 

24.3531 

10.094 

203.4 

21.5145 

10.032 

197.7 

18.2143 

9.957 

188.2 

16.6258 

9.919 

188.2 

14.8533 

9.875 

181.8 

14.6597 

9.870 

184.7 

13.4571 

9.839 

181.5 

11.4848 

9.787 

175.0 

7.6669 

9.679 

166.0 

4.1001 

9.568 

155.9 

.8064 

9.453 

148.2 

.4145 

9.439 

147.1 

T = 250.00  K 


31.6827 

10.970 

377.8 

27.9692 

10.917 

3u4.3 

24.5552 

10.866 

349.7 

21.0180 

10.811 

338.7 

19.2675 

10.784 

327.8 

17.5178 

10.755 

327.7 

T26 


Table  III  3 (Continued) 


P,  hPa 

p, 

/LiPa* 

15.5723 

10.723 

314.7 

14.0358 

10.697 

316.8 

11.6838 

10.656 

302.7 

10.7614 

10.640 

304.7 

7.3892 

10.578 

295.0 

3.1925 

10.496 

279.9 

T = 200.00  K 


31.2686 

11.687 

748.8 

28.1272 

11.654 

731.7 

24.8552 

11.619 

711.1 

21.3494 

11.580 

682.6 

17.6364 

11.538 

658.1 

16.8536 

11.528 

643.6 

14.0934 

11.496 

638.0 

12.6896 

11.479 

609.6 

10.7472 

11.455 

620.2 

7.3582 

11.412 

601.4 

3.0687 

11.356 

583.8 

T27 


Table  III  3 (Continued) 


P,  MPa  p,  iiiol'L”^  ?j,|jPa*s 

T = 180.00  K 


31.4359 

11.978 

1093.6 

28.0430 

11.947 

1064.2 

24,4960 

11.913 

1017.0 

21.0666 

11.880 

979.9 

17.6586 

11.846 

945.5 

13.9976 

11.808 

909.3 

10.5922 

11.772 

884.3 

7.4037 

11.737 

858.0 

3.3993 

11.691 

836.7 

T = 160.00  K 


34.1087 

12.291 

1842.7 

30.6682 

12.263 

1802.2 

27.1278 

12.235 

1688.6 

25.6699 

12.223 

1694.5 

21.9056 

12.191 

1656.0 

18.3818 

12.160 

1584.2 

14.8377 

12.128 

1505.8 

11.1760 

12.095 

1469.4 

7.9163 

12.064 

1441.0 

4.4249 

12.029 

1381.8 

1.6752 

12.002 

1346.0 

T28 


Table  III  3 (Continued) 


P,  MPa  p,  mol'L"^  ^^^.yjpa^s 

T = 150.00  K 


32.3674 

12.423 

2459.5 

29.3326 

12.400 

2412.6 

25.9258 

12.374 

2341.6 

22.6434 

12.348 

2246.8 

18.5270 

12.315 

2145.4 

14.9878 

12.286 

2077.1 

14.8903 

12.285 

2116.8 

11.5634 

12.256 

1972.2 

10.8120 

12.250 

1981.1 

7.7792 

12.223 

1906.1 

7.3933 

12.219 

1944.2 

3.3096 

12.182 

1890.5 

T = 140.00  K 

34.4848 

12.584 

3540.8 

31.4535 

12.564 

3482.9 

27.5545 

12.536 

3226.3 

24.0140 

12.510 

3156.6 

20.3420 

12.483 

3014.1 

16.9883 

12.457 

2957.9 

15.1308 

12.443 

2966.6 

12.4894 

12.422 

2826.7 

9.7485 

12.400 

2722.5 

T29 


Table  III  3 


(Continued) 


P,  MPa 

p,  inol'L  ^ 

yi*  uPa. 

7.5270 

12.381 

2640.4 

4.9825 

12.359 

2597.0 

2.2351 

12.335 

2512.6 

T 

= 135.00 

K 

32.8042 

12.647 

4404.4 

29.3833 

12.624 

4310.2 

26.0042 

12.600 

4143.9 

22.4957 

12.576 

3994.0 

18.9878 

12.550 

3878.8 

15.4475 

12.524 

3671.2 

11.9075 

12.496 

3523.9 

7.2201 

12.459 

3352.1 

2.9658 

12.424 

3137.1 

T 

= 130.00 

K 

32.3245 

12.717 

5358.4 

29.1835 

12.697 

5239,7 

25.8556 

12.675 

5049.6 

22.4324 

12.651 

4849.0 

19.0082 

12.628 

4708.6 

15.4097 

12.602 

4515.2 

11.8763 

12.576 

4359.5 

7.7479 

12.544 

4214.7 

2.7303 

12.504 

4028.7 

T30 


Table  III  3 (Continued) 


P.  MPa 

p,  mol*L 

t2» 

T = 125.00  K 

33.4321 

12.798 

7088.0 

30.2429 

12.778 

6803.0 

27.1840 

12.759 

6562.1 

23.9961 

12.738 

6244.6 

20.3856 

12.714 

6078.9 

16.9765 

12.691 

5742.6 

13.1952 

12.665 

5553.6 

9.8415 

12.640 

5366.3 

6.4430 

12.615 

5138.3 

2.8559 

12.587 

4800.9 

T = 120.00  K 

33.1409 

12.870 

9302.9 

30.2300 

12.853 

8772.1 

27.0404 

12.834 

8604.2 

23.9612 

12.815 

8269.6 

20.2683 

12.791 

7925.3 

16.8658 

12.769 

7516.3 

13.1879 

12.744 

7371.2 

10.4375 

12.725 

7198.2 

7.5484 

12.704 

6908.2 

3.2824 

12.673 

6471.3 

T31 


Table  III  ^ 


1 


ities  of  saturated  liquid  isobutane 


T,  K 

-1 

P,  Irtol'L 

uPd»S 

115.00 

12.733 

8188.7 

115.00 

12.733 

8224.1 

116.00 

12.716 

7652.4 

118.00 

12.683 

6919.0 

120.00 

12.647 

6334.7 

125.00 

12.567 

4694.8 

130.00 

12.481 

3751.5 

135.00 

12.401 

2999.1 

140.00 

12.315 

2447.3 

140.00 

12.315 

2412.9 

145.00 

12.235 

2091.6 

150.00 

12.150 

1747.5 

155.00 

12.070 

1480.5 

160.00 

11.984 

1272.1 

160.00 

11.984 

1287.2 

160.00 

11.984 

1301.5 

170.00 

11.817 

975.3 

180.00 

11.651 

829.3 

180.00 

11.651 

803.5 

180.00 

11.651 

797.5 

190.00 

11.483 

651. 9 

200.00 

11.313 

539.0 

200.00 

11.313 

531.5 

T32 
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Table  III  4 (Continued) 


T,  K 

7^,  uPa.s 

220.00 

10.969 

396.2 

220.00 

10.969 

393.1 

240.00 

10.614 

301.8 

240.00 

10.614 

298.7 

250.00 

10.431 

262.2 

260.00 

10.243 

231.9 

260.00 

10.243 

234.2 

280.00 

9.853 

183.9 

280.00 

9.853 

184.9 
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Fig,  I 1 Arraingement  of  visual  cell,  pressure  transducer, 

thermostated  bath  and  thermometer  for  the  oberserva- 
tion  of  meniscus  disappearance  in  pure  fluids  and 
mixtures. 


^4 


vent 


Fig.  I 2 Arrangement  of  visual  cell,  thermal  compressor  (T-C), 

injector  pump  (I-P)  and  pressure  transducer  (P)  used  for 
density  and  composition  determinations  in  binary  fluid 
mixtures.  During  oberservations,  the  entire  manifold 
is  filled  with  mercury,  and  the  cell  in  part. 


Fig.  II  1 The  critical  line  of  isobutane-isopentane 

mixtures  plotted  against  the  mole  fraction  x of 
isopentane:  (a)  critical  temperature,  (b)  critical 
density  and  (c)  critical  pressures.  The  points 
represent  the  experimental  data  of  Table  I 5 and  the 
curves  are  calculated  from  the  Leung-Grif f iths  equations 
(Sec.  3 of  Part  II). 


(a)  (b) 


Fig.  II  2 Corapariaon  between  experimental  and  calculated  dew- 

bubble  curve  data  at  various  mole  fractions  x of  isopentane 
(a)  coexistence  curves  near  the  top  of  the  dome  in  the 
temperature  density  plane,  and  (b)  vapor- pres sure 
curves  in  the  pressure  - temperature  plane.  The  points 
represent  the  experimental  data  of  Table  I 7 and  the  curves 
the  values  calculated  from  the  Leung-Griff iths 
equations  (See.  3 of  Part  II).  - 
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Fig.  II  3a  Comparison  between  experimental  dew-bubble  curve  data, 

the  classical  equation  of  state,  and  the  Holdover-Rainwater 
model,  in  the  temperature  - (mass)  density  plane.  The 
dashed  lines  are  the  critical  lines,  the  narrow  solid 
lines  are  the  pure  fluid  coexistence  curves,  and  the  bold 
solid  lines  are  calculated  from  the  Holdover-Rainwater 
model.  The  points  near  critical  represent  the  data  of 
Table  I 7,  and  chose  away  from  critical  represent  the 
predictions  of  the  classical  equation  of  state.  Compositions 
are:  0.2  0;  0.35  A;  0.5079  □;  0.65  7 (Sec.  5 of  Parc  II). 


Figo  II  3b  Comparison  between  experimental  dew"bubMe  curve  dataj,  the 
classical  equation  of  state*  and  the  Moldover-Rainwater 
model,  pressure“temperature  plane.  The  dashed  lines  are  the 
critical  lines,  the  narrow  solid  lines  are  the  pure  fluid 
coexistence  curves,  and  the  bold  solid  lines  are  calculated 
from  the  Moldover-Rainwater  model.  The  points  near  critical 
represent  the  data  of  Table  I 7,  and  chose  away  from  critical 
represent  the  predictions  of  the  classical  equation  of  state. 
Compositions  are:  0.20;  0.35  A;  0.5079  □;  0.65  V (Sec,  5 of 
Part  II) . 
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Fig.  Ill  1 Viscosity  of  compressed  gaseous  and  liquid  methane  as 
a function  of  density  along  various  isotherms. 


13.000 


Fig.  Ill  2 Viscosity  of  compressed  gaseous  and  liquid  ethane  as 
a function  of  density,  (p^  is  the  critical  density; 

is  the  density  of  the  liquid  at  the  triple  point 
temperature, ) 
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Fig.  Ill  3 Viscosity  of  compressed  gaseous  and  liquid  propane 
a function  of  density  cilong  various  isotherms. 
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Fig.  Ill  M Viscosity  of  compressed  gaseous  and  liquid  50.115  mol 
nitrogen  + 49.885  mol  % methane  as  a function  of  densi 
along  various  isotherms. 
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Fig,  III  5 Viscosity  of  compressed  gaseous  auid  liquid  0.50217  mol  % 
methcuie  + 49.783  mol  % ethane  as  a function  of  density. 
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Fig.  Ill  6 Viscosity  of  compressed  gaseous  and  liquid  normal 
butane  at  saturation  aind  along  various  isotherms. 


VISCOSITY,  ixPa  s 


DENSITY,  mol 


Fig.  Ill  7 Viscosity  of  compressed  gaseous  and  liquid  isobutauie 
at  saturation  and  along  various  isotherms. 
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Fig.  Ill  8 Comparison  of  the  viscosity  of  methane  at  100  K with 
the  correlation  of  Ref.  51. 
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Fig.  Ill  9 Comparison  of  the  viscosity  of  ethane  with  that  pre- 
dicted by  the  generalized  corresponding-states  model. 
0,  290  K;  250  K;  □,  200  K;  A,  150  K;  ^,130  K;  7 , 
120  K; f , 110  K;  O,  100  , saturated  liquid. 
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Fig.  Ill  10  Comparison  of  the  viscosity  of  propane  with  that 
predicted  by  the  generalized  corresponding-states 
model.  0,  300  K;  140  K;  □,  110  K;  A,  100  K; 

95  K;  , 90  K.  •,  saturated  liquid. 
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Fig.  Ill  11  Comparison  of  the  viscosity  of  a 50.115  mol  % 

nitrogen — 49.885  mol  % methane  mixture  with  that 
predicted  by  the  generalized  corresponding-states 
model.  0,  155  K;  140  K;  □,  120  K;  A,  110  K; 

V , 100  K. 
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Fig.  ill  12  Viscosity  of  coinprossed  liquid  nitrogsn  + mothane 

mixtures  and  their  pure  components  as  a function  of 
composition  at  fixed  molar  density  (28.0  mol  •liter”"' ) 
and  at  fixed  temperature  (100  K).  The  present  mea- 
surements, ; mole  fraction  average  of  the  pure 
component  viscosities,  — — — ; generalized  corre- 
sponding states  model, . 


Fig.  Ill  13  Comparison  of  the  viscosity  of  compressed  liquid 
methane  and  ethane  with  the  generalized  corre- 
sponding-states model. 
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Fig.  Ill  14  Comparison  of  the  viscosity  of  a 35.528  mol  % 

methane-65.472  mol  % ethane  mixture  with  the  gener- 
alized corresponding-states  model. 
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Fig.  Ill  15  Dependences  of  the  viscosities  of  saturated  liquid 
normal  butane  and  saturated  liquid  isobutane  on 
density.  Normal  butane;  ©,  measured  viscosities; 
generalized  corresponding  states  model.  Isobutane 
□ , measured  viscosities; — generalized  corre- 

sponding states  model. 
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Fig.  Ill  16  Differences  between  measured  normal  butane  visco- 
sities, measured  isobutane  viscosities  and  a gener- 
alized corresponding  states  model. 

•,300  K;  •,  250  K;4,  200  K;A,  1 80  K;  ^,160  K;  □, 

150  K;  4 ,140  K;  V , 135  K;  A,  130  K;  V ,125  K;  & 
,120  K;  o,  saturated  liquid. 


(VISCOSITY)"’ (/igcm-’-s‘’r’ 


^ III  17  Dependence  of  the  fluidity  of  ethane  on  molar  volume 
and  temperature. 
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Fig.  Ill  18 


Dependence  of  the  fluidity  of  propane  on  molar 
volume  and  temperature. 
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Fig.  Ill  19  Dependence  of  the  fluidities  of  saturated  and 

compressed  liquid  normal  butane  and  isobutane  on 
molar  volume  and  temperature. 
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butane  from  a Hildebrand- type  correlating  equation. 
•,  300  K;  a,  250  K;4,  200  K;  A,  1 80  K;  f,  160  K;  □, 
150  K;  , 140  K.  o,  saturated  liquid. 
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Fig.  Ill  21  Percentage  departures  of  the.  viscosity  of  isobutane 
I from  a Hildebrand- type  correlating  equation.  • , 

I 300  K;  a,  250  K;^,  200  K;  4,  1 80  K;  i 60  K;  □, 

j 150  K;  > , 140  K;  V , 135  K;  A,  130  K;  r , 125  K; 

) . 3 » 120  K.  o,  saturated  liquid. 
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